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MULTI-DIMENSIONAL JOINT SEARCHER AND 
CHANNEL ESTIMATORS 

BACKGROUND 



[0001] This application is related to the following United States Patent applications, all 

5 simultaneously filed herewith: United States Patent Application Serial 10/ , 

(attorney docket 2380-795), entitled "Spatial Joint Searcher And Channel Estimators"; 

United States Patent Application Serial 10/ , (attorney docket 2380-796), entitled 

"Temporal Joint Searcher And Channel Estimators", United States Patent Application 

Serial 10/ , (attorney docket 2380-797), entitled "Spatio-Temporal Joint Searcher 

10 And Channel Estimators", all of which are incorporated by reference herein. 

[0002] FIELD OF THE INVENTION 

[0003] The present invention pertains to wireless telecommunications, and particularly 
to apparatus and method for determining a channel estimate for use in reconstructing 
data symbols transmitted over a channel. 

15 [0004] RELATED ART AND OTHER CONSIDERATIONS 

[0005] A wireless telecommunications unit typically includes both a transmitter and 
receiver for communicating with other wireless telecommunications units over a 
communication link. For wireless communications, the communication link is typically 
over an air interface (e.g., radio frequency interface). As used herein, a "wireless 

20 telecommunications unit" with its "wireless telecommunications receiver" can be 
included in a network node (e.g., a radio access network node such as a base station 
node, also called Node-B) or a terminal. Such "terminals" include mobile terminals 
such as user equipment units (UEs), which have also been called mobile stations, and 
include by way of example mobile telephones ("cellular" telephones), laptops with 

25 mobile termination. Thus, terminals can be, for example, portable, pocket, hand-held, 
computer-included, or car-mounted mobile devices which communicate voice and/or 
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data with radio access network. Alternatively, the terminals can be fixed wireless 
devices, e.g., fixed cellular devices/terminals which are part of a wireless local loop or 
the like. 

[0006] As shown simply in Fig. 32, a wireless telecommunications system includes a 
5 transmitting antenna 2300T and a receiving antenna 2300R. Channel 2302 describes 
the relation between the transmitting antenna 2300T and the receiving antenna 2300R, 
including the wireless interface. A signal, typically modulated into pulses, is 
transmitted over channel 2302 from transmitting antenna 2300T to receiving antenna 
2300R. The signal can comprise a "symbol" or a string of series of symbols, depicted 

10 as "m" in Fig. 32. The signal can carry user data and/or certain control data (e.g., a 
pilot bit or pilot sequence). The signal m as transmitted by the transmitting antenna 
2300T is convoluted with a channel impulse response h of the channel, so that the 
received signal at the receiving antenna 2300R is m * h (e.g., m convoluted with h). 
The received signal m*his applied to base band processing functionality 2304 of the 

15 receiver where the received signal undergoes radio frequency processing. The data 
portions of the received signal are applied to a detector 2306, which may be, for 
example, a demodulator such as a RAKE receiver. 

[0007] Most modern detectors attempt to recover a symbol estimate m from the 
received signal m * h . To do so, most sophisticated detectors expect to receive a 
20 "channel estimate" for use in modeling the channel over which the signal was 
transmitted. The accuracy of this channel estimate influences the accuracy and 
performance of the detector in estimating the actual symbol received over the channel. 

[0008] The modeling of the channel (which is necessary for most detectors) is 
facilitated by the control data, often in the form of a pilot bit or pilot sequence, which is 

25 transmitted by the transmitter. The control data, hereinafter referenced as "pilot data" 
for simplicity, is of a known or recognizable format or pattern. The pilot data is 
typically transmitted periodically by the transmitter source, and thus receipt of 
repetitions of the pilot data can be expected at the receiver at successive intervals. In 
view of factors such as relative motion of the transmitter and receiver, the successive 

30 intervals are not necessarily constant. The pilot data can be transmitted simultaneously 
with, or otherwise interspersed with, the user data. 
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[0009] In order to utilize the pilot data, wireless receivers typically include both a 
searcher and a channel estimator, such as searcher 2308 and channel estimator 2310 
shown in Fig. 32. For control data, the received signal m* his applied to the searcher 
2308, which determines a time of arrival (TOA). The time of arrival is then applied to 
5 the channel estimator 2310, which uses the time of arrival to determine the channel 
estimate h and then provides the channel estimate h to the detector 2306. Using the 
channel estimate h , the detector develops its estimate of the symbol, e.g., m . 

[00010] The receiver may receive the an original signal (e.g., short pulse signal) from 
the transmitter source through open space over a single, direct propagation path. 
10 Alternatively, in another environment having obstacles or other surfaces, the receiver 
may receive the same original signal over multiple propagation paths. In the multiple 
path case, the received signal appears at the receiver as a stream of pulses, each pulse 
having a different time delay in view of the corresponding propagation multipath over 
which the signal travelled, as well as possibly different amplitude and phase. 

15 [0001 1] Multipaths are created in a mobile radio channel by reflection of the signal 

from obstacles in the environment such as buildings, trees, cars, people, etc. Moreover, 
the mobile radio channel is dynamic in the sense that it is time varying because of 
relative motion affecting structures that create the multipaths, or due to movement of 
structures and objects in the surroundings (even if the transmitter and receiver are 

20 fixed). For a signal transmitted over a time varying multipath channel, the received 
corresponding multiple paths vary in time, location, attenuation, and phase. 

[00012] Some wireless telecommunications receivers capitalize upon the existence of 
the multipaths in order to achieve various advantages. Such receivers typically operate 
on the baseband signal to search for and identify the strongest multipaths along with 

25 their corresponding time delays. The receiver has a filter which operates on a power 
delay profile of the signal. The power delay profile can be conceptualized as a time- 
averaged refinement or other derivation of the channel impulse response. The searcher 
attempts to locate peaks in the power delay profile, each peak corresponding to arrival 
of a wavefront of the signal from a respective multipath. In many searchers the peaks 

30 also correspond to a channel tap of the filter. 
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[00013] A channel estimate ft as applied to the detector therefore comprises a set of 
both time of arrivals (TO A) and complex channel coefficients, each pair of TOA and 
channel coefficients being associated with one of the arriving wavefronts. In other 
words, each arriving wavefront has a pair of members in the set, e.g., a TOA and a 
5 channel coefficient. The channel coefficients thus actually form a channel impulse 
response vector, so that the terms "channel coefficient" and "channel coefficients" as 
used hereinafter should be understood to refer to a channel impulse response vector. If 
there is only one wavefront, there is only one TOA and one channel coefficient in the 
set (one channel coefficient in the channel impulse response vector) . But for plural 
10 arriving wavefronts, there are a corresponding plurality of TO As and channel 

coefficients. Ideally, the channel estimate h should provide as good an estimate of the 
channel impulse response as possible, thereby increasing performance of the detector as 
the detector makes its estimate m of the transmitted symbol m. 

[00014] The channel estimate is then supplied to the detector, such as RAKE type of 
15 demodulators. A RAKE demodulator typically allocates a number of parallel 

demodulators (called RAKE fingers) to the strongest multipath components of the 
received multipath signal as determined by the multipath search processor. In a 
wideband code division multiple access (WCDMA) radio access network, the outputs 
of each of the RAKE fingers are diversity-combined after corresponding delay 
20 compensation to generate a "best" demodulated signal that considerably improves the 
quality and reliability of the radio communications system. 

[00015] Conventionally, wireless telecommunications receivers first use their searchers 
to ascertain time of arrival of a wavefront. Subsequently, after the time of arrival has 
been determined by the searcher, the channel estimator utilizes the time of arrival to 
25 calculate a channel coefficient, which expresses both amplitude and phase of the signal. 

[00016] Some wireless telecommunications units have more than one antenna for 
receiving a same signal. In the prior art, the searcher attempts to locate peaks in the 
power delay profile for each antenna separately. In other words, for each antenna the 
searcher works more or less independently. See, for example, US Patent Publication 
30 US 2002/0048306, which is incorporated herein by reference. As such, the prior art 
searchers are essentially one dimensional. 
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[00017] As indicated above, the performance of a wireless receiver is considerably 
dependent upon the accuracy of the peak determination, i.e., time of arrival 
determination, performed by the searcher. The better the peak determination of the 
searcher, the better will be the overall performance of the receiver (e.g., less error rate). 

5 But in many instances it may be difficult for a searcher to find an actual peak in a 
power delay profile. As mentioned previously, in many searcher algorithms the peak 
corresponds to a channel tap. With such difficulty there is considerable risk of 
incorrectly choosing a peak. Moreover, it can then be difficult to estimate the actual 
channel tap value. Channels with low signal to noise ratios (SINRS) are particularly 

10 susceptible to these difficulties. 

[00018] What is needed, therefore, and an object of the present invention, is provision 
of apparatus and method for providing an improved channel estimate for a wireless 
telecommunications receiver. 



BRIEF SUMMARY 

15 [00019] A wireless communication receiver comprises an antenna structure and a joint 
searcher and channel estimator. The antenna structure acquires dimensionally 
differentiated signals which are concurrently utilized by the joint searcher and channel 
estimator for determining both a time of arrival and channel coefficient. The wireless 
communication receiver can be either a mobile terminal or a network node (e.g., a radio 

20 access network node such as a base station node, also called Node-B). 

[00020] In some embodiments, the antenna structure comprises an array of plural 
antennas. The signals acquired by different antennas of the array are dimensionally 
differentiated with regard to a spatial dimension. As such, the joint searcher and 
channel estimator of these embodiments are "spatial" joint searcher and channel 
25 estimators. The time of arrival and the channel coefficient are essentially concurrently 
determined by the joint searcher and channel estimator. The time channel coefficient is 
a composite channel coefficient which takes into consideration channel impulse 
responses for channels associated with each of the plural antennas in the antenna array. 



[00021] In other embodiments, the antenna structure comprises an antenna which 
30 provides signals for each of successive sets of pilot data received by the antenna as the 
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dimensionally differentiated signals. The signals acquired by the antenna are 
dimensionally differentiated with regard to a temporal dimension. As such, the joint 
searcher and channel estimator of this embodiment is a ''temporal" joint searcher and 
channel estimator. 

5 [00022] In one of its aspects, the joint searcher and channel estimator comprises an 
antenna signal matrix; a matrix analyzer; and a channel estimate generator. Complex 
values indicative of the dimensionally differentiated signal received in a sampling 
window are stored in the antenna signal matrix as a function of a sampling window 
time index and a dimensional differentiation index. In differing embodiments, the 
10 matrix analyzer can employ non-parametric or parametric techniques. 

[00023] In non-parametric embodiments, the matrix analyzer can take the form a 
correlator which performs to produce a vector or matrix. The correlator locates value(s) 
in the antenna signal matrix for use in determining the time of arrival and the channel 
coefficient. Specifically, the correlator locates 0 "tones" in the antenna signal matrix 
15 corresponding to arriving wavefronts. The channel estimate generator uses the value(s) 
located by the correlator to generate the time of arrival and the channel coefficient. 

[00024] In locating values in the antenna signal matrix corresponding to an arriving 
wavefront, the correlator considers a dimensional reception vector formed from the 
antenna signal matrix with respect to a sampling window time index. The dimensional 

20 receptivity vector has a frequency related to a difference between phase components of 
complex values of the dimensional receptivity vector, there being plural possible 
frequencies for the dimensional receptivity, the plural possible frequencies being 
represented by a frequency index. In the spatial embodiments, the plural possible 
frequencies correspond to plural possible directions of arrival (DOA) of a wavefront. 

25 In the temporal embodiments, the plural possible frequencies correspond to plural 

possible frequency shifts for an arriving wavefront. The frequency shifts are primarily 
attributable to doppler shifts, but can also include frequency shifts in the transmitter and 
receiver oscillators. For simplification, such frequency shifts are hereinafter referred to 
as "doppler shifts" or "doppler frequency shifts". 

30 [00025] In one example non-parametric type correlator implementation for a joint 

searcher and channel estimator, for each combination of plural possible frequencies and 
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plural time indexes the correlator calculates a Fast Fourier Transform according to the 
expression Y(n,t) = FFT(n,X(:,t)). In this expression, t is the sampling window time 
index; X(:,t) is the complex antenna matrix, with : representing all dimensional 
differentiation indexes for one sampling window time index; and n is the frequency 
5 index. For a CDMA wireless communication system, the correlator calculates the 
expression Y(n,t) = X C j *FFT(n,X(:,t)), j = 1,K, wherein Cj is a coding sequence 
symbol value j, and K is a length of the coding sequence. In spatial embodiments, the 
dimensional differentiation indexes are the antenna indexes, and each of the plural 
possible frequencies for the dimensional receptivity vector represents a different 
10 possible direction of arrival of the arriving wavefront. In temporal embodiments, the 
dimensional differentiation indexes are pilot set indexes, and each of the plural possible 
frequencies for the dimensional receptivity vector represents a different possible 
doppler shift values for an arriving wavefront. 

[00026] The non-parametric type correlator output comprises Y(n,t). The analyzer 
15 determines a maximum absolute value |Y(n,t)| max . The analyzer uses a sampling 
window time index t_max at which |Y(n,t)| max occurs as the time of arrival of the 
arriving wavefront. In the spatial embodiments, the analyzer uses the a frequency index 
n_max at which |Y(n,t)| max occurs as the direction of arrival of the arriving wavefront 
and obtains an amplitude for the arriving wavefront by dividing |Y(n,t)| max by a number 
20 of antennas comprising the antenna array. In the spatial embodiments, the direction of 
arrival is used to derive the channel coefficient. In the temporal embodiments, the 
analyzer uses the a frequency index n_max at which |Y(n,t)| max occurs as the doppler 
shift for the arriving wavefront and obtains an amplitude for the arriving wavefront by 
dividing |Y(n,t)| max by a number of sets of pilot data in the series. In the temporal 
25 embodiments, the doppler shift is used to derive the channel coefficient. 

[00027] In parametric embodiments, the matrix analyzer can take the form of a 
estimator which estimates a parameter in the signal, and thereby generates a parameter 
output estimation value (e.g., a parametic estimation output vector). The parametic 
estimation output vector typically includes spatial frequency, spatial amplitude, 
30 temporal frequency and temporal amplitude or spatio-temporal frequency and spatio- 
temporal amplitude. The parametic estimation output vector could include other 
parameters. The parametric matrix analyzer generates the parametric estimation output 
vector. The channel estimate generator uses the parametric estimate output vector to 
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generate the time of arrival and the channel coefficient. The channel estimate generator 
uses amplitude and frequency values of elements of the parametric estimation output 
vector to determine the time of arrival and channel coefficient of an arriving wavefront. 

[00028] For the spatial embodiments, the parametric estimation output vector 
5 generated by the parametric estimator has a sampling window time index and a spatial 
parameter value for each time index. The spatial parameter includes spatial frequency 
and spatial amplitude. For an element of the parametric estimation output vector 
having a sufficiently high absolute value, i.e. the spatial amplitude, the channel estimate 
generator uses a sampling window time index for an element of the parametric 
10 estimation output vector having a sufficiently high spatial amplitude value to determine 
the time of arrival of the arriving wavefront. The direction of arrival of the arriving 
wavefront is enhanced from the spatial frequency parameter value of the identified time 
of arrival. 

[00029] For the temporal embodiments, the parametric estimation output vector 
15 generated by the parametric estimator has a sampling window time index and a 
temporal parameter value for each time index. The temporal parameters includes 
temporal frequency and temporal amplitude. For an element of the parametric 
estimation output vector having a sufficiently high absolute value, the channel estimate 
generator uses a sampling window time index for an element of the parametric 
20 estimation output vector having a sufficiently high absolute value, i.e. the temporal 

amplitude, to determine the time of arrival of the arriving wavefront. The doppler shift 
of the arriving wavefront is the temporal frequency parameter value of the identified 
time of arrival. 

BRIEF DESCRIPTION OF THE DRAWINGS 

25 [00030] The foregoing and other objects, features, and advantages of the invention will 
be apparent from the following more particular description of preferred embodiments as 
illustrated in the accompanying drawings in which reference characters refer to the 
same parts throughout the various views. The drawings are not necessarily to scale, 
emphasis instead being placed upon illustrating the principles of the invention. 
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[00031] Fig. 1 is a schematic view of an example, generic wireless telecommunications 
receiver which includes a joint searcher and channel estimator. 

[00032] Fig. 2 A and Fig. 2B are schematic views of differing example embodiments of 
spatial joint searcher and channel estimators, each shown with an antenna array. 

5 [00033] Fig. 3 is a diagrammatic view illustrating a signal emanating from a 

transmitting antenna along three separate multipaths to an antenna array of a wireless 
telecommunications receiver. 

[00034] Fig. 4 is a diagrammatic view of a wavefront travelling toward an antenna 
array. 

10 [00035] Fig. 5A and Fig. 5B are diagrammatic views depicting signals obtained upon 
arrival of a wavefront at an antenna array. 

[00036] Fig. 6 is a diagrammatic view of an antenna signal matrix. 

[00037] Fig. 7 is a flowcharting showing representative basic steps performed by a 
matrix analyzer and channel estimate generator of an example embodiment of a spatial 
15 joint searcher and channel estimator, with the matrix analyzer using a non-parametric 
analysis technique. 

[00038] Fig. 8A, Fig. 8B, Fig. 8C(1), Fig. 8C(2), and Fig. 8C(3) and are diagrammatic 
views depicting results of a comparative operational evaluation contrasting 
performance of a spatial joint searcher and channel estimator with a conventional 
20 searcher. 

[00039] Fig. 9A is a diagrammatic view of an antenna signal matrix; an antenna weight 
vector; and a non-parametric output estimation vector. 

[00040] Fig. 9B is a diagrammatic view of an antenna signal matrix and a parametric 
output estimation vector. 

25 [00041] Fig. 10 is a flowcharting showing representative basic steps performed by 

matrix analyzer and channel estimate generator of an example embodiment of a spatial 
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joint searcher and channel estimator, with the matrix analyzer using a parametric 
analysis technique. 

[00042] Fig. 1 1 is a diagrammatic view illustrating a coherent combination of signal 
outputs by a joint searcher and channel estimator. 

5 [00043] Fig. 12 is a schematic view for illustrating how an antenna weight vector 
facilitates the coherent combination illustrated in Fig. 1 1 . 

[00044] Fig. 13A is a schematic view of an example embodiment of a temporal joint 
searcher and channel estimator shown with an antenna array, the temporal joint 
searcher and channel estimator comprising a matrix analyzer which employs a non- 
10 parametric analysis technique. 

[00045] Fig. 13B is a schematic view of an example embodiment of a temporal joint 
searcher and channel estimator shown with an antenna array, the temporal joint 
searcher and channel estimator comprising a matrix analyzer which employs a 
parametric analysis technique. 

15 [00046] Fig. 14 is a diagrammatic view depicting a sequence of sets of pilot data and 
user data received by a receiver which utilizes a temporal joint searcher and channel 
estimator, as well as an antenna signal matrix utilized by the temporal joint searcher 
and channel estimator. 

[00047] Fig. 15 is a flowcharting showing representative basic steps performed by a 
20 matrix analyzer and channel estimate generator of an example embodiment of a 

temporal joint searcher and channel estimator, with the matrix analyzer using a non- 
parametric analysis technique. 

[00048] Fig. 16A is a diagrammatic view of an antenna signal matrix; a doppler weight 
vector; and a non-parametric estimation output vector for a temporal joint searcher and 
25 channel estimator. 



[00049] Fig. 16B is a diagrammatic view of an antenna signal matrix and a parametric 
estimation output vector for a temporal joint searcher and channel estimator. 
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[00050] Fig. 17 is a flowcharting showing representative basic steps performed by a 
matrix analyzer and channel estimate generator of an example embodiment of a 
temporal joint searcher and channel estimator, with the matrix analyzer using a 
parametric analysis technique. 

5 [00051] Fig. 18A is a schematic view of an example embodiment of a spatio-temporal 
joint searcher and channel estimator shown with an antenna array, the spatio-temporal 
joint searcher and channel estimator comprising a matrix analyzer which employs a 
non-parametric analysis technique. 

[00052] Fig. 18B is a schematic view of an example embodiment of a spatio-temporal 
10 joint searcher and channel estimator shown with an antenna array, the spatio-temporal 
joint searcher and channel estimator comprising a matrix analyzer which employs a 
parametric analysis technique. 

[00053] Fig. 19 is a diagrammatic view depicting a sequence of sets of pilot data and 
user data received by a receiver which utilizes a combined spatial/temporal joint 
15 searcher and channel estimator, as well as an antenna signal matrix utilized thereby. 

[00054] Fig. 20 is a flowcharting showing representative basic steps performed by a 
matrix analyzer and channel estimate generator of an example embodiment of a spatio- 
temporal joint searcher and channel estimator, with the matrix analyzer using a non- 
parametric analysis technique. 

20 [00055] Fig. 21 is a diagrammatic view of an antenna signal matrix; a doppler weight 
and antenna weight vector; and a non-parametric estimation output vector for an 
example embodiment of a spatio-temporal joint searcher and channel estimator which 
operates in a three dimensional essentially concurrent mode. 

[00056] Fig. 22A and Fig. 22B are diagrammatic views depicting operation of a first 
25 alternative implementation of a non-parametric, sequential spatio-temporal joint 
searcher and channel estimator. 

[00057] Fig. 23 describes the procedure of non-parmetric approach for spatio-temporal 
sequenced method where the spatial processing is followed by the temporal processing. 
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[00058] Fig. 24A and Fig. 24 Bare diagrammatic views depicting operation of a second 
alternative implementation of a non-parametric, sequential spatio-temporal joint 
searcher and channel estimator. 

[00059] Fig. 25 describes the procedure of non-parmetric approach for spatio-temporal 
sequenced method where the temporal processing is followed by the spatial processing. 

[00060] Fig. 26 is a diagrammatic view of an antenna signal matrix and a parametric 
estimation output vector for an example embodiment of a spatio-temporal joint searcher 
and channel estimator. 

[00061] Fig. 27 is a flowcharting showing representative basic steps performed by a 
matrix analyzer and channel estimate generator of an example embodiment of a spatio- 
temporal joint searcher and channel estimator, with the matrix analyzer using a 
parametric analysis technique. 

[00062] Fig. 28A and Fig. 28B are diagrammatic views depicting operation of a first 
alternative implementation of a parametric, sequential spatio-temporal joint searcher 
and channel estimator. 

[00063] Fig. 29 describes the procedure of the parametric approach for spatio-temporal 
sequenced method where the spatial processing is followed by the temporal processing. 

[00064] Fig. 30A and Fig. 30B are diagrammatic views depicting operation of a second 
alternative implementation of a parametric, sequential spatio-temporal joint searcher 
and channel estimator. 

[00065] Fig. 31 describes the procedure of the parametric approach for spatio-temporal 
sequenced method where the temporal processing is followed by the spatial processing. 

[00066] Fig. 32 is a schematic view of a conventional wireless telecommunications 
receiver. 

DETAILED DESCRIPTION OF THE DRAWINGS 

[00067] In the following description, for purposes of explanation and not limitation, 
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specific details are set forth such as particular architectures, interfaces, techniques, etc. 
in order to provide a thorough understanding of the present invention. However, it will 
be apparent to those skilled in the art that the present invention may be practiced in 
other embodiments that depart from these specific details. In other instances, detailed 
5 descriptions of well-known devices, circuits, and methods are omitted so as not to 
obscure the description of the present invention with unnecessary detail. Moreover, 
individual function blocks are shown in some of the figures. 

[00068] Fig. 1 shows an example, generic wireless telecommunications receiver 20 
which, as mentioned before, can be included in a network node or a terminal, e.g. 
10 mobile terminal. The wireless telecommunications receiver 20 includes an antenna 
structure or array 22; a joint searcher and channel estimator 24; a detector 26; and, a 
timing and control unit 28. Optionally, as depicted by broken line, the receiver 20 may 
include a code sequence generator 30. 

[00069] As broadly employed herein, the antenna array 22 can comprise one or more 
15 antenna elements. Signal(s) from the antenna array 22 are applied both to joint searcher 
and channel estimator 24 and detector 26. The signal(s) from the antenna array 22 
comprise a channel impulse response vector if the antenna array 22 comprises more 
than one antenna element. 

[00070] In the likely event that the signal(s) have been encoded by, e.g., a spreading 
20 code or the like, both joint searcher and channel estimator 24 and detector 26 are 

connected to operate in conjunction with code sequence generator 30. The timing and 
control unit 28 generates timing (e.g., synchronization) and control signals which are 
provided to detector 26 and to joint searcher and channel estimator 24. 

[00071] It will be appreciated that the receiver may include, e.g., downstream from the 
25 antenna array, certain radio frequency processing functionality and radio frequency 
demodulating functionality, so that the signals applied to joint searcher and channel 
estimator 24 and detector 26 are baseband signals. The illustrated structure of wireless 
telecommunications receiver 20 of Fig. 1 thus essentially concerns processing of the 
baseband signal(s). 

30 [00072] Various non-limiting, representative examples of differing embodiments of 
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joint searcher and channel estimators are described below. Ensuing descriptions of 
operation of wireless telecommunications receivers with these differing embodiments 
are premised on certain assumptions. Some of these assumptions are related to a 
channel model which conceptualizes electromagnetic fields as arriving in a discrete 
5 number of wavefronts at the wireless telecommunications receiver, and particularly 
arriving at one or more antenna elements which may be employed in the antenna array 
22. 

[00073] As used herein, a "sampling window" comprises consecutive time slots (or, in 
a CDMA system, for example, "chips") obtained from a given antenna and analyzed by 

10 a joint searcher and channel estimator. As described in more detail hereinafter, 

embodiments of joint searcher and channel estimators operate upon an antenna signal 
matrix formed from plural sampling windows. In some embodiments, hereinafter 
referenced as "spatial" joint searcher and channel estimators, the antenna signal matrix 
is formed from sampling windows obtained from plural antennas. In other 

15 embodiments, hereinafter referenced as "temporal" joint searcher and channel 

estimators, the antenna signal matrix is formed with respect to a single antenna, but 
formed from sampling windows obtained by that antenna for successive sets of pilot 
data (occurring over time). In yet other embodiments, hereinafter referenced as spatio- 
temporal joint searcher and channel estimators, the antenna signal matrix is formed 

20 both spatially and temporally. 

[00074] For purposes of the technology described herein, the antenna array 22 is 
conceptualized as acquiring "dimensionally differentiated" signals. The joint searcher 
and channel estimator essentially concurrently uses the dimensionally differentiated 
signals provided by the antenna array for determining, for each arriving wavefront, both 

25 a time of arrival (TO A) and a channel coefficient. For the spatial joint searcher and 
channel estimator, wherein the antenna structure comprises an array of plural antennas 
having spaced apart or spatially separated antenna elements, the signals acquired by 
different antennas of the array are dimensionally differentiated with regard to a spatial 
dimension. For the temporal joint searcher and channel estimator, wherein the antenna 

30 structure comprises an antenna which provides signals for each of successive sets of 
pilot data received at separated time intervals, the signals acquired by the antenna are 
dimensionally differentiated with regard to a temporal or time dimension. For the 
spatio-temporal joint searcher and channel estimator, having both the antenna structure 
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comprising an array of plural antennas and one or more antennas receiving the 
successive sets of pilot data, the signals acquired by the antenna are dimensionally 
differentiated with regard both to a spatial dimension and a temporal or time dimension. 

[00075] The joint searcher and channel estimators are, in some instances, said to 
5 perform a "concurrent" determination of time of arrival and some other quantity, e.g., 
direction of arrival or doppler shift frequency. In this sense "concurrent" means that the 
quantities or determinations could be derived in parallel from a result of an outcome- 
determinative operation, e.g. a non-parametric technique such as a Fast Fourier 
Transform or a parametric technique. 

10 [00076] SPATIAL JOINT SEARCHERS/ESTIMATORS 

[00077] In some embodiments, the joint searcher and channel estimator 
contemporaneously processes the signals from plural antennas over a sampling window 
in order to determine both time of arrival and the channel coefficient. In these 
embodiments the joint searcher and channel estimator is essentially a two dimensional 

15 unit, with a second dimension being a spatial dimension imparted by the spacing of the 
plural antennas of the array. This spatial dimension, which involves essentially 
simultaneous and concurrent processing of signals from the plural antennas for the 
array in order to determine the time of arrival and channel coefficient, bestows on these 
embodiments of the joint searcher and channel estimator the distinction of being a 

20 "spatial" joint searcher and channel estimator. 

[00078] The spatial joint searcher and channel estimator may take differing 
embodiments and have differing implementations. In one example, illustrative 
embodiment the joint searcher and channel estimator includes a non-parametric type 
correlator (e.g., a correlator which performs a Fast Fourier Transform (FFT) 
25 calculation). In another example, illustrative embodiment the joint searcher and 
channel estimator utilizes a parametric approach 

[00079] Fig. 2A illustrates one example embodiment of a spatial joint searcher and 
channel estimator 24-2A which uses a non-parametric technique for determining time 
of arrival and channel estimate, as well as an associated example antenna array 22-2 A. 
30 The antenna array 22-2A includes, by way of non-limiting example, four antenna 
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elements 22-2A-1 through 22-2A-4. While the antenna elements 22-2A-1 through 22- 
2A-4 are shown as forming a uniform linear array (ULA), it should be understood that 
antenna configurations other than a uniform linear are possible, and that the number of 
antenna elements in the antenna array may vary (e.g., the number of antenna elements 
5 is not limited to four). 

[00080] There are coherency requirements for the antenna elements of antenna array 
22-2A, and for antenna elements for all other plural antenna arrays described herein. 
The coherency requirement can be fulfilled by the plural antenna elements being 
synchronized. Alternatively, even if the plural antenna elements are not synchronized, 
10 but their phase differences are known, the coherency requirement can be fulfilled by 
compensating for the known phase difference. 

[00081] The complex baseband signals obtained from the antenna elements are each 
applied to joint searcher and channel estimator 24-2A, as well as to a detector (not 
illustrated in Fig. 2A). The joint searcher and channel estimator 24-2A comprises an 
antenna signal matrix handling unit 40-2 A. In one particular example manifestation, 
antenna signal matrix handling unit 40-2A includes antenna signal matrix generator 42- 
2A and antenna signal matrix memory 44-2A. A matrix analyzer, which for the non- 
parametric technique of Fig. 2A can be a correlator 50-2A, operates on complex values 
stored in antenna signal matrix memory 44-2 A. The correlator 50-2 A preferably 
comprises a filter. The correlator 50-2A generates certain output values, which may be 
stored, e.g., in correlator output value memory 52-2 A. The joint searcher and channel 
estimator 24-2 A further comprises a channel estimate (CE) generator 60-2 A. In the 
illustrated example embodiment, the channel estimate (CE) generator 60-2A comprises 
a correlator output analyzer 62-2 A and a detector interface 64-2 A. The detector 
interface 64-2A generates, for each wavefront, both a time of arrival (TOA) and a 
channel coefficient (CC). In Fig. 2A, the time of arrival and channel coefficient output 
by detector interface 64 are applied to the detector on lines 66-2A and 68-2A, 
respectively. 

[00082] In Fig. 2A, and other embodiments described herein, the transmitted 
30 electromagnetic signal is assumed to arrive at the receiver in a number of discrete 
electromagnetic wavefronts. A number of discrete electromagnetic wavefronts is 
presumed in order to accommodate the multipath phenomena discussed above. For 
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example, Fig. 3 illustrates a signal emanating from a transmitting antenna 70 along 
three separate multipaths P u P 2 , and P 3 to antenna array 22. Each multipath has its 
individual amplitude, and accordingly has an associated complex number "a" of the 
baseband signal and a time delay x. For example, multipath Pi has associated complex 
5 number a x and associated time delay x^ multipath P 2 has associated complex number a 2 
and associated time delay x 2 ; and so forth. As illustrated in Fig. 3, multipath Pi is a 
relatively direct path between transmitting antenna 70 and antenna array 22; while 
multipath P 2 and multipath P 3 are reflected off obstacles 72 2 and 72 3 , respectively. 
Thus, the time delay Xi for multipath Pi is shorter than the time delay x 2 for multipath 
10 P 2 , which in turn is shorter than the time delay x 3 for multipath P 3 . Similarly, barring 
other phenomena, it would be expected that the complex number ai for multipath Pi is 
greater than the complex number a 2 for multipath P 2 , and so forth. 

[00083] For sake of discussion, the electromagnetic wavefronts are assumed to be plane 
("planar") electromagnetic wave fronts, such as the single wavefront 76 illustrated in 
15 Fig. 4 as traveling toward the antenna array. In all embodiments described herein it 
should be understood that the wavefronts need not be planar wavefronts, but that any 
other known form of wavefront may be considered in similar manner. Moreover, it 
should be kept in mind that Fig. 4 represents arrival of only one wavefront, but 
typically plural wavefronts are incident on an antenna array. 

20 [00084] As further shown by Fig. 4, due to incidence of an individual wavefront, the 
output (e.g., signal) from each antenna element has its version of the complex number 
for the wavefront. For example, for a wavefront for the first multipath Pi of Fig. 3, the 
antenna element 22-1 outputs a complex number ai_i, antenna element 22-2 outputs a 
complex number ai_ 2 , and so forth. The numbers are complex, and in the particular case 

25 that (1) the antenna elements are identical; (2) there is coherency, and (3) the plane 
wave has constant amplitude within the width of array, the absolute values of the 
numbers are the same. Furthermore, with respect to the same arriving wavefront, each 
antenna detects the arriving signal as having a phase. For example, for the wavefront 
for the first multipath Pi of Fig. 3, the output of antenna element 22-1 has a phase 0 M , 

30 the output of antenna element 22-2 has a phase 0i_ 2 , and so forth. 

[00085] Signals obtained upon arrival of a wavefront at a uniform linear array (ULA) 
antenna array are illustrated both in Fig. 5 A and in Fig. 5B. Fig. 5 A particularly shows, 
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for each of four antennas 22-1 through 22-4, plane wave propagation over the antenna 
elements for a fixed time (chip) index, and resulting respective output pulses 78 (e.g., 
output pulses 78i through 78 4 ). For each corresponding antenna, Fig. 5B shows the 
pulse as a complex number and with the argument,0, of the complex number. The 
5 argument (6) corresponds to the phase of the received signal. The rate at which the 0 
values change (e.g., the rate at which the phase rotates) over time is known as the phase 
rotation speed, or frequency. The phase rotation for the wavefront with this array of 
antennas is depicted by the increasing angular value of 0 through the range of 01, 0 2 , 03, 
0 4 , and thus the frequency is the rate of change of this angular value over time. The 
10 phase rotation speed is constant. The speed of the linear phase propagation is 
dependent on the direction of arrival (DOA) of the incident wavefront. 

[00086] In the joint searcher and channel estimator 24-2A of Fig. 2A, the antenna 
matrix handling unit 40-2A samples the complex baseband signals from each antenna 
element. Using the sampled complex baseband signals, antenna signal matrix generator 
15 42-2A generates an antenna signal matrix such as antenna signal matrix 80 illustrated in 
Fig. 6. The antenna signal matrix 80 may be stored in any convenient fashion, such as 
antenna matrix memory 44-2 A. 

[00087] The antenna signal matrix 80 is a two dimensional functionally dependent 
matrix. In other words, complex samples are stored in antenna signal matrix 80 as a 

20 function of two different indexes. For the antenna signal matrix 80 shown in Fig. 6, a 
first index is a sampling window time index, illustrated along the X axis of Fig. 6. For 
embodiments which utilize spreading codes or similar codes, the first index may be, for 
example, a chip index. Thus, the sampling window time index points to a time in the 
sampling window relative to a start of the sampling window. In the antenna signal 

25 matrix 80 of Fig. 6, a second index, shown along the Y axis, is an antenna index (which 
serves as a dimensional differentiation index). The antenna index points to a different 
row of the antenna signal matrix 80, each row being associated with a different antenna 
element in antenna array 22. Fig. 6 shows four rows in antenna signal matrix 80 for 
consistency with the previous examples of an antenna array comprising four antenna 

30 elements. It is reiterated, however, that the number of antennas in an antenna array, and 
thus the number of rows in antenna signal matrix 80 and the maximum value of the 
antenna index, can vary from receiver to receiver, and that the choice of four antenna is 
only illustrative for sake of example. 
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[00088] The antenna signal matrix 80 is conceptualized as storing "dimensionally 
differentiated" signals acquired from the antenna array. For the spatial joint searcher 
and channel estimator, wherein the antenna structure comprises an array of plural 
antennas having spaced apart or spatially separated antenna elements, the signals 
5 acquired by different antennas of the array are dimensionally differentiated with regard 
to a spatial dimension. That is, for a given column of antenna signal matrix 80, the 
values in each row are dimensionally differentiated in the sense that they are acquired 
from different antenna elements which are separated in a spatial dimension in view of 
the separate physical placement of each antenna element with respect to other antenna 
10 elements of the array. 

[00089] For sake of simplicity, the complex values stored in antenna signal matrix 80, 
including the complex values obtained from the antennas, are not illustrated in Fig. 6. 
Such complex values would be illustrated in a third dimension, e.g., out of the plane of 
Fig. 6. The antenna signal matrix 80 includes both complex white noise and (for the 
15 sake of the present illustration) a complex sample for at least one wavefront (planar or 
other known shape). As stored in antenna signal matrix 80, the wavefronts have known 
phase (temporal, non-coherent detection), and are modulated code sequences. 

[00090] In conjunction with the antenna signal matrix 80 of Fig. 6, and particularly a 
WCDMA case in which spacing of the antenna elements in the antenna array is not too 
20 far apart, the plane wavefront arriving at the antenna array can be considered to arrive 
in the same sampling window time index (or chip index). 

[00091] The complex values stored for each column of the antenna signal matrix 80 of 
Fig. 6 can be conceptualized as a dimensional receptivity vector. That is, a dimensional 
receptivity vector is formed with respect to a single sampling window time instance and 

25 with complex values from each of the plural antennas of the antenna array. Each 

element taken from a unique row of antenna signal matrix 80 has a different phase in 
the manner of the differing 0 values illustrated in Fig. 5. As received by the differing 
antenna elements, for the spatial joint searcher and channel estimator the change in 
phase over time is the frequency for the dimensional receptivity vector. If the wave 

30 arrives e.g. straight ahead, the angles could be the same. The phase rotation speed, or 
frequency, of the dimensional receptivity vector, for the sampling window time 
instance can be interpreted as the direction of arrival (DO A). Thus, each dimensional 



receptivity vector corresponds to a separate direction of arrival. There are plural 
possible frequencies for the dimensional receptivity vector, each of the plural possible 
frequencies corresponding to a different possible direction of arrival (DOA) of a 
wavefront. For the non-parametric techniques herein employed, the plural possible 
5 frequencies can be a continuous range of frequencies. For sake of differentiating the 
plural possible frequencies, the plural possible frequencies are each represented by a 
frequency index. 

[00092] The channel estimate generator 60-2A (see Fig. 2A) seeks to develop a 
"composite" channel estimate based on the complex values stored in antenna signal 

10 matrix 80. At this point it should be appreciated that, since the antenna array 22-2A has 
plural antenna elements, there are a corresponding plurality of channels through which 
wavefronts are received, and accordingly there could also be a separate channel impulse 
response or separate channel estimate for each of the plural channels. But by storing 
the complex samples in antenna signal matrix 80 in the manner aforedescribed, and by 

15 concurrently finding the time of arrival (TO A) and channel coefficients over the entire 
antenna signal matrix 80, the channel estimate generator 60-2A provides a channel 
estimate which encompasses all channels for all antenna elements and for this reason is 
known as a "composite" channel estimate. 

[00093] The composite channel estimate comprises, as mentioned before, a time of 
20 arrival (TO A) and channel coefficient for each arriving wavefront in the sampling 

window (e.g., a channel coefficient mapped to a time of arrival (TOA)). Therefore, the 
channel estimate may comprise a set (of one or more) pairs of data, each pair including 
a time of arrival (TOA) and channel coefficient. The task for correlator 50-2A is thus 
to locate a value or "tone" in antenna signal matrix 80 that best corresponds to an 
25 arriving wavefront, e.g., to locate a value or tone for each arriving wavefront in the 
sampling window. 

[00094] The task of locating a value or "tone" in antenna signal matrix 80 that best 
corresponds to an arriving wavefront can be accomplished by various techniques, 
including both parametric and non-parametric techniques. A Fast Fourier Transform 
30 (FFT) technique as discussed below is just one representative and illustrative example 
non-parametric type of correlator which can be utilized. 
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[00095] Fig. 7 depicts example basic steps performed by an example correlator 50-2A 
and correlator output analyzer 62-2A in conjunction with the Fast Fourier Transform 
(FFT) calculation. As step 7-1, the correlator 50-2A of Fig. 2A calculates Expression 1. 

Y(n,t) = FFT(n,X(:,t)) Expression 1 

5 In Expression 1, t is the sampling window time index; X(:,t) is the complex antenna 
matrix (with the colon ":" representing all antenna indexes for one sampling window 
time index); and, n is the frequency index. Each FFT calculation is thus a one 
dimensional FFT calculation on the baseband signal, and corresponds to a specific 
direction of arrival (as depicted by the frequency index) and set of antenna weights 
10 which, in practice, are the FFT weights. 

[00096] The output of correlator 50-2A, i.e., the Y(n,t) values computed using 
Expression 1, are stored as correlator output values. The correlator output values can 
be stored, for example, in the correlator output value memory 52-2 A of Fig. 2A. 

[00097] The correlator output analyzer 62-2A of channel estimate (CE) generator 60- 
15 2A searches the correlator output values and (as step 7-2) determines therefrom a 
maximum absolute value |Y(n,t)| max . This maximum absolute value |Y(n,t)| max is 
utilized by correlator output analyzer 62-2A to determine both the direction of arrival 
(DOA) and time of arrival (TO A) for an arriving wavefront seen in the sampling 
window. In particular, as step 7-3 correlator output analyzer 62-2 A chooses a sampling 
20 window time index t__max at which |Y(n,t)| max occurs to be the time of arrival of the 
arriving wavefront. In addition, as step 7-4 correlator output analyzer 62-2A chooses 
the frequency index n_max at which |Y(n,t)| max occurs to represent the direction of 
arrival (DOA) of the arriving wavefront. The frequency index corresponds to a 
direction of arrival (e.g., 0). An amplitude for the arriving wavefront is determined as 
25 correlator output analyzer 62-2A divides |Y(n,t)| max by the number of antennas 
comprising the antenna array (as step 7-5). 

[00098] Expression 1 and the steps of Fig. 7 represent a generic non-parametric FFT 
calculation. In a CDMA-specific situation which utilizes a coding generator (such as 
coding generator 30 of Fig. 1), a comparable FFT calculation can be made using a 
30 refinement of Expression 1 which appears as Expression 2. 
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Y(n,t) = S Cj *FFT(n,X(:,t)), j = 1,K Expression 2 

Expression 2 is understood from Expression 1, it being further mentioned that Cj is a 
coding sequence symbol value j; and K is a length of the coding sequence. 

[00099] As a result of operation of joint searcher and channel estimator 24-2 A, an 
5 accurate channel estimate can be provided to the detector as a spatial signature. The 
spatial signature includes the time of arrival (TOA), as well as the direction of arrival 
(DOA) and amplitude. As explained below, the channel coefficient (CC) for each 
wavefront is derived from the direction of arrival (DOA) and amplitude. The time of 
arrival (TO A) and channel coefficient (CC) are applied to the detector as represented by 
10 lines 66-2A and 68-2A, respectively, in Fig. 2A. 

[000100] As mentioned above, the channel coefficient (CC) for each wavefront is 
derived from the direction of arrival (DOA) and amplitude. Recall that at step 7-4 
correlator output analyzer 62-2 A chose the frequency index n_max at which |Y(n,t)| max 
occurs to represent the direction of arrival (DOA) of the arriving wavefront, with the 
15 chosen frequency index corresponding to a direction of arrival (e.g., 0). The channel 
impulse response vector (i.e., array propagation vector) x is therefore generated by 
detector interface 64-2A in accordance with Expression 3 (for identical isotropic 
antenna elements). 

[000101] x =[l,e akd * sin9) , e akd ^ sin9) , . . . e akd * (K1)sine) ]*C Expression 3 

20 [000102] In Expression 3, j is the conventional imaginary notation; k = 2*ttA,; d is a 
spacing distance between elements of the antenna array; A, is the wavelength of the 
received/transmitted electro-magnetic signal: (f*A =c) and, K is the antenna element 
index (illustrated as antenna numbers Al, A2, A3, A4 in FIG. 9 A, for example) In 
Expression 3, C is a complex constant in which |C| = |FFTjiiax|/number of antennas; 

25 the argument of C, i.e., arg(C) = arg(FFTjnax), wherein |FFT_max is the FFT value 
computed at step 7-1 of Fig. 7. 

[000103] In the foregoing description it is the role of channel estimate (CE) 
generator 60-2A, and particularly detector interface 64-2A, to generate both a time of 
arrival (TOA) and a channel coefficient (CC), the channel coefficient being derived 
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from the direction of arrival, e.g., as above described in conjunction with Expression 3. 
In an alternate implementation of this and other embodiments described herein, the 
detector itself (such as detector 26 illustrated in Fig. 1), upon receiving the time of 
arrival (TOA) and direction of arrival (DOA) for each arriving wavefront, may have the 
5 intelligence to compute the channel coefficient for each wavefront from the 

corresponding direction of arrival (DOA) information. In such case, the time of arrival 
and direction of arrival are output by detector interface 64 to the detector. 

[000104] Thus, considering the aspects above discussed, the joint searcher and 
channel estimator 24-2A looks in a discrete number of possible directions of arrival, 

10 and picks the direction of arrival with the highest correlation (highest absolute value). 
A comparative operational evaluation was preformed to illustrate the efficacy of a joint 
searcher and channel estimator such as joint searcher and channel estimator of 24-2A of 
Fig. 2A. A first scenario of the comparative operational evaluation involved a 
conventional searcher functioned essentially in prior art fashion for a sampling window. 

15 In so doing, with respect to each antenna for the sampling window the conventional 
search merely picked the time (e.g., chip) which had the greatest absolute value. In 
other words, the signals from each antenna were processed separately. A second 
scenario of the comparative operational evaluation was performed in the manner above 
described with respect to the joint searcher and channel estimator 24-2A of Fig. 2 and 

20 Expression 1 . The same signal was applied in both scenarios to an antenna array 

having eight antenna elements. The length of the sampling window for both scenarios 
was twenty chips, and a coding sequence of { 1 } was utilized (e.g., only one of the chips 
contained the signal, the remainder of the chips contained complex white noise). 

[000105] Fig. 8 A illustrates the first scenario which utilized the conventional 
25 searcher. In contrast, Fig. 8B illustrates the spatial joint searcher and channel estimator 
24-2A of Fig. 2A utilized for the second scenario. The superiority of the second 
scenario (and thus the spatial joint searcher and channel estimator) is evident by a 
comparison of Fig. 8 A and Fig. 8B, due to the higher SNR for the signal of interest in 
Fig. 8B . In the second scenario, it is much easier to pick out the tone or value for the 
30 arriving wavefront. For the second scenario, Fig. 8C(1) shows the absolute value of the 
complex channel impulse response taps; Fig. 8C(2) shows phase errors of the complex 
channel impulse response taps; and Fig. 8C(3) shows the detected time of arrival. 
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[000106] Whereas the joint searcher and channel estimator of Fig. 2A includes a non- 
parametric type matrix analyzer, e.g., a correlator (e.g., a filter which performs a Fast 
Fourier Transform (FFT) calculation), in other example embodiments the matrix 
analyzer of the joint searcher and channel estimator implements parametric techniques. 
5 As does the Fig. 2A embodiment, the spatial joint searcher and channel estimator 24-2B 
of Fig. 2B (which uses a parametric technique) is shown along with its associated 
example antenna array 22-2B. Again by way of example, antenna array 22-2B includes 
four antenna elements 22-2B-1 through 22-2B-4. The signals obtained from the 
antenna elements are each applied to joint searcher and channel estimator 24-2B, as 
10 well as to a detector (not illustrated in Fig. 2B). 

[000107] Similar to the earlier described embodiment, joint searcher and channel 
estimator 24-2B can comprise an antenna signal matrix handling unit 40-2B, which in 
turn comprises antenna signal matrix generator 42-2B and antenna signal matrix 
memory 44-2B, which function much in the manner previously described. For 
15 example, the complex baseband values stored in antenna signal matrix memory 44-2B 
can also be conceptualized as matrix 80, and as such has a sampling window time 
index. The antenna signal matrix 80 has been previously discussed in conjunction with 
Fig. 6, and is now also discussed with reference to Fig. 9A for sake of expounding the 
joint searcher and channel estimator24-2B of Fig. 2B. 

20 [000108] The joint searcher and channel estimator 24-2B further comprises a matrix 
analyzer, e.g., parametric estimator 51 -2B, which utilizes a parametric technique. In 
addition, in similar manner as the preceding embodiment, joint searcher and channel 
estimator 24-2B comprises a channel estimate generator 60-2B which has parametric 
estimation output vector analyzer 62-2B and a demodulator interface 64-2B. Basic 

25 steps performed by parametric estimator 51-2B and parametric estimation output vector 
analyzer 62-2B of the joint searcher and channel estimator 24-2B of Fig. 2B are 
illustrated in Fig. 10. 

[000109] For each sampling window time index of the antenna signal matrix 80, as step 
10-1 the parametric estimator 51-2B estimates, e.g. two parameters at each time instant: 
30 a spatial frequency parameter parameter and a spatial amplitude parameter. The spatial 
frequency parameter estimates the frequency the incident waves creates when arriving 
at the ULA. The spatial amplitude parameter estimates the amplitude of this frequency. 
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The spatial frequency parameter and spatial amplitude parameter are considered to be a 
parameter pair and in Fig. 9B, they are illustrated as one parameter per sample along 
the sampling time index. The parameters can be calculated by an appropriate strategy 
or goal criteria, e.g., by a minimum mean square error technique (MMSE). 

5 [0001 10] As step 10-2, parametric estimation output vector analyzer 62-2B finds certain 
"qualifying" values in parametric estimation output vector, i.e. high or maximum values 
of the spatial amplitude parameter. The qualifying values can be, for example, values 
whose absolute values are sufficiently high or are a maximum. Each qualifying value 
of parametric estimation output vector 90 can correspond to an arriving wavefront for 

10 the sampling window. 

[0001 11] For each qualifying value, as step 10-3 the parametric output estimation 
vector analyzer 62-2B chooses a time of arrival (TO A) as corresponding to the 
sampling window time index t for the qualifying value, e.g., the time index at which the 
maximum/qualifying absolute value of the parametric estimation output vector occurs. 

15 [0001 12] Similarly, for each qualifying value, as step 10-4 the analyzer 62-2B chooses 
a direction of arrival (DOA) as the spatial frequency parameter value at the time of 
arrival, decided in 10-3. 

[0001 13] As step 10-5, the parametric estimation output vector analyzer 62-2B 
determines the amplitude as the value of the spatial amplitude value divided by the 
20 number of antenna elements in the array. 

[0001 14] The joint searcher and channel estimator 24-2B thus looks for an 
optimum direction, and prepares a channel estimate which can be provided to the 
detector as a spatial signature. The spatial signature includes the direction of arrival 
(DOA) and amplitude. The channel coefficient (CC) for each wavefront is derived 
25 from the direction of arrival (DOA) and amplitude in the manner explained above with 
reference to Expression 3. The time of arrival (TOA) and channel coefficient (CC) are 
applied to the detector as represented by lines 66-2B and 68-2B, respectively, in Fig. 
2B. 

[0001 15] It should be understood from the foregoing that information indicative of 
30 more than one incident wavefront may be seen in a sampling window. For example, 
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with reference to the parametric estimation output vector 90 of Fig. 9B, the parametric 
estimation output vector analyzer 62-2B may see other high numbers and for each of 
those high numbers which qualify, an arriving wavefront may be ascertained. For 
example, if there were two high numbers, then the channel impulse response may 
5 reflect two arriving wavefronts. For each of the two arriving wavefronts the joint 
searcher and channel estimator would pick out both a time of arrival (TOA) and 
direction of arrival (DOA), as well as amplitude, which are mapped to two different 
channel coefficients, with these two different channel coefficients forming part of the 
channel estimate. 

10 [0001 16] Fig. 4 showed that a wavefront individually reached each of four example 
antenna elements of an antenna array, providing a different antenna output (complex 
baseband signal) for each antenna element. For example, the output of antenna 
element 22-1 with the complex vector a^ (and phase 0i_i); the output of antenna 
element 22-2 is the complex vector a U2 (and phase 0m), and so forth. The linear 

15 combination, of the complex antenna baseband signal and the antenna weight vectors 
Wj have the effect of a summation, or coherent combination in the time and space 
domain, shown as summation function 100 in Fig. 12. 

[0001 17] The coherent combination facilitated by the antenna weight vectors Wj is 
illustrated in Fig. 11. In the example case of the four antenna elements shown in Fig. 

20 12, the effect of weight belonging to antenna index 2, here denoted asW 2 is to rotate the 
output of antenna element 22-2 so that its phase 9i_ 2 lines up with the phase 9 m of the 
output of antenna element 22-1, in the manner shown in Fig. 1 1 . Similarly, the effect of 
weight W 3 is to rotate the output of antenna element 22-3 so that its phase 0!_ 3 lines up 
with the phase 0 m of the output of antenna element 22-1. The effect of weight W 4 is 

25 to rotate the output of antenna element 22-4 so that its phase 0i_ 4 lines up with the phase 
6m of the output of antenna element 22-1 . For simplicity, Fig. 1 1 ignores noise 
considerations, which tend to make the resultant vector less than straight. Note that in 
the preciding paragraph, the weight vectors are denoted with Wi, where i denotes the 
antenna index of the weight vector W, which is denoted without index. 

30 [0001 18] In the spatial joint searcher and channel estimators, the SINR for finding 
channel taps (peaks) should be proportional to the number of antenna elements 
comprising the array. The operation of the spatial joint searcher and channel estimators 
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can be adapted to take into consideration channel variations over time, e.g., spatial 
variations in the environment (e.g., in the sending and receiving antennas). 

[0001 19] The non-parametric FFT-type correlator and the parametric techniques 
illustrated above, e.g., by Fig. 2A and Fig. 2B, respectively, are only two example 
5 techniques for finding the values or "tones" in antenna signal matrix 80 which are 

associated with arriving wavefronts. Other parametric approaches are described by or 
understood from Stocia, Petre and Moses, Randolph, Introduction To Spectral Analysis, 
ISBN-013-258419-0, Prentice Hall, which is incorporated by reference in its entirety, 
particularly Chapter 4 thereof. 

10 [000120] The spatial joint searcher and channel estimator and techniques of 

operation thereof as described above are suitable for any receiver unit which has plural 
receiving antennas. Thus, the spatial joint searcher and channel estimator is particularly 
well suited for, but not limited to, a base station which has plural antennas. Also 
encompasses are mobile terminals which have plural antennas. 

15 [000121] TEMPORAL JOINT SEARCHERS/ESTIMATORS 

[000122] In other embodiments, the joint searcher and channel estimator 
contemporaneously processes the signals received at an antenna element from plural, 
successive sets of pilot data (each set of pilot data being received in its own sampling 
window) in order to determine both time of arrival and the channel coefficient. In so 

20 doing, the joint searcher and channel estimator takes into consideration a doppler shift 
or frequency shift (the terms "doppler shift" and "frequency shift" being used 
interchangeably in conjunction with the description of the temporal joint searcher and 
channel estimator). The frequency shift is primarily attributable to a doppler shift, but 
can also include a frequency shift in the transmitter and receiver oscillators. For 

25 simplification, such frequency shifts are hereinafter referred to as "doppler shifts" or 
"doppler frequency shifts". 

[000123] The doppler shift can be occasioned by movement such as relative 
movement of one of the transmitter and the receiver (for example, by movement of a 
mobile terminal), or movement of a signal path-affecting object or structure in the 
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surroundings (which can cause a doppler shift even for a fixed transmitter and fixed 
receiver). 

[000124] In providing the channel estimate, the joint searcher and channel 
estimator essentially concurrently considers plural signals (e.g., plural sets of pilot data) 
5 received by the antenna element. The joint searcher and channel estimator applies the 
channel coefficient and the time of arrival to a detector which provides, e.g., a symbol 
estimate. 

[000125] In these embodiments, the joint searcher and channel estimator is 
essentially a two dimensional unit, with a second dimension being a temporal 

10 dimension imparted by the time intervals at which the successive sets of pilot data 
arrive. This temporal dimension, which involves essentially simultaneous and 
concurrent processing together of signals received at the antenna element from each of 
the plural sets of pilot data, bestows on these embodiments of the joint searcher and 
channel estimator the distinction of being a "temporal" joint searcher and channel 

15 estimator. 

[000126] The temporal joint searcher and channel estimator may take differing 
embodiments and have differing implementations. In one example, illustrative 
embodiment the temporal joint searcher and channel estimator includes a non- 
parametric type correlator (e.g., a correlator which performs a Fast Fourier Transform 
20 (FFT) calculation). In another example, illustrative embodiment the temporal joint 
searcher and channel estimator utilizes a parametric approach. 

[000127] Fig. 13A illustrates one example embodiment of a spatial joint searcher and 
channel estimator 24-13A which uses a non-parametric technique for determining time 
of arrival and channel estimate, as well as an associated example antenna array 22- 13 A. 
25 In the example of Fig. 13 A, the antenna array 22-13A is shown as having one antenna 
element 22-13A-L As explained hereinafter, complex baseband signals obtained from 
the same antenna element (e.g., antenna element 22-13A-1) upon receipt of each of 
successive sets of pilot data (as hereinafter described) are each applied to joint searcher 
and channel estimator 24- 13 A, as well as to a detector (not illustrated in Fig. 13 A). 

30 [000128] The joint searcher and channel estimator 24-1 3 A comprises an antenna signal 
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matrix handling unit 40- 13 A. In one particular example manifestation, antenna signal 
matrix handling unit 40-1 3 A includes antenna signal matrix generator 42-13A and 
antenna signal matrix memory 44- 13 A. A matrix analyzer, which for the non- 
parametric technique of Fig. 2A can be correlator 50-1 3 A, operates on complex values 

5 stored in antenna signal matrix memory 44- 13 A. The correlator 50-1 3 A preferably 
comprises a filter. The correlator 50-13A generates certain output values, which may 
be stored, e.g., in correlator output value memory 52- 13 A. The joint searcher and 
channel estimator 24-13A further comprises a channel estimate (CE) generator 60- 13 A. 
In the illustrated example embodiment, the channel estimate (CE) generator 60-1 3 A 

10 comprises a correlator output analyzer 62-1 3 A and a detector interface 64- 13 A. The 
detector interface 64-13A generates, for each wavefront, a channel estimate which 
includes both a time of arrival (TO A) and a channel coefficient (CC). In Fig. 13 A, the 
time of arrival and channel coefficient output by detector interface 64-13A are applied 
to the detector on lines 66- 13 A and 68- 13 A, respectively. 

15 [000129] As shown in Fig. 14, the temporal joint searcher and channel estimators such 
as joint searcher and channel estimator 24-13A of Fig. 13A watch the channel response 
from an antenna (e.g., antenna 22-13-1) for sets of pilot data which are interpersed or 
otherwise transmitted with other data (e.g., user data). For sake of simplicity, it is 
assumed that each set of pilot data is received in a separate sampling window. Such 

20 need not be the case, however, as differing sets of pilot data can be received 

simultaneously if the different streams are, e.g., code multiplexed. Merely as an 
illustrative example, Fig. 14 shows four sets of pilot data, i.e., pilot sets Tl - T4, 
interspersed with user data and received at unique global times (as indicated by the "T" 
axis in Fig. 4). 

25 [000130] Each set of pilot data is typically in a different frame from another set of pilot 
data. For example, pilot set Tl may be in frame 1 ; pilot set T2 may be in frame 1 1 ; 
pilot set T3 may be in frame 21; etc. "Frame transmission interval 1 ' refers to the time 
between two successive frames which contain pilot data. The time between two 
successive frames which contain pilot data is typically specified by a standard or other 

30 specification. 

[000131] Fig. 14 thus reflects the typical periodic transmission of pilot data by the 
transmitter source, and also the expected receipt of repetitions of the pilot data at the 
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receiver at successive intervals. In view of factors such as relative motion of the 
transmitter and receiver, the successive intervals between differing sets of pilot data are 
not necessarily constant. 

[000132] As further shown in Fig. 14, an antenna matrix handling unit (such as antenna 
5 matrix handling unit 40-13A of the Fig. 13A embodiment) samples the signals received 
by the antenna element for each of the successive sets of pilot data, i.e., for pilot sets Tl 
- T4. Using the sampled signals, antenna signal matrix generator 42-13A generates an 
antenna signal matrix such as antenna signal matrix 110 illustrated in Fig. 14. The 
antenna signal matrix 110 may be stored in any convenient fashion, such as antenna 
10 matrix memory 44- 1 3 A. 

[000133] The antenna signal matrix 1 10 is a two dimensional functionally dependent 
matrix. In other words, complex samples are stored in antenna signal matrix 1 10 as a 
function of two different indexes. For the antenna signal matrix 110 shown in Fig. 14, 
a first index is a sampling window time index, illustrated along the X axis of Fig. 14. 

15 For embodiments which utilize spreading codes or similar codes, the first index may be, 
for example, a chip index. Thus, the sampling window time index points to a time in 
the sampling window relative to a start of the respective sampling window. In the 
antenna signal matrix 1 10 of Fig. 14, a second index, shown along the Y axis, is a pilot 
set index (which serves as a dimensional differentiation index). The pilot set index 

20 indicates which one of the sets of pilot data the sample was obtained. In other words, a 
pilot set index = Tl indicates that the sample was obtained from pilot set Tl ; a pilot set 
index = T2 indicates that the sample was obtained from pilot set T2; and so on as 
depicted by the arrows which connect the matrix 110 with the received signal depiction 
with its illustrative successive sets of pilot data. As can be seen, the pilot set index 

25 points to a different row of the antenna signal matrix 110, each row being associated 
with a different set of pilot data. 

[000134] Fig. 14 shows four rows in antenna signal matrix 1 10 for consistency with the 
illustrated example wherein the antenna signal matrix encompasses four successive sets 
of pilot data. The number of sets of pilot data subsumed in a given antenna signal 
30 matrix, and thus the maximum value of the pilot set index, can vary from receiver to 
receiver, so that the present example's choice of four sets of pilot data is only 
illustrative for sake of example. In general, the choice of the number of sets of pilot 
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data to be apprehended simultaneously by a temporal joint searcher and channel 
estimator depends on how quickly the doppler is expected to change. The number of 
taps/incident waves depends on the multipath. In other words in an open space we have 
one direct path and thus only one channel/tap coefficient in the channel impulse 
5 response. 

[000135] The antenna signal matrix 1 10 is also conceptualized as storing 
"dimensionally differentiated" signals acquired from a single antenna element of the 
antenna array. For the temporal joint searcher and channel estimator, wherein the 
antenna structure comprises an antenna which provides signals for each of successive 
10 sets of pilot data received at separated time intervals, the signals acquired by the 

antenna are dimensionally differentiated with regard to a temporal or time dimension. 
For example, the signals acquired by the antenna are dimensionally differentiated by 
being acquired in differing frame transmission intervals. 

[000136] For sake of simplicity, the complex values stored in antenna signal matrix 
15 110, including the complex values obtained from the antennas, are not illustrated in Fig. 
14. Such complex values would be illustrated in a third dimension, e.g., out of the 
plane of Fig. 14. The antenna signal matrix 110 includes both complex white noise and 
(for the sake of the present illustration) a complex sample for at least one wavefront 
(planar or other known shape). The wavefronts have known phase (temporal, non- 
20 coherent detection), and are modulated code sequences. 

[000137] The complex values stored for each column of the antenna signal matrix 
1 10 of Fig. 14 can be conceptualized as a dimensional receptivity vector. That is, a 
dimensional receptivity vector is formed with the complex values taken with respect to 
a same single sampling window time index for each of the sets of pilot signals included 

25 in the sampling window (e.g., for sets Tl - T4 in Fig. 14). Each element taken from a 
unique row of antenna signal matrix 110 has a different phase in the manner of the 
differing 0 values illustrated in Fig. 5. As received by the differing antenna elements, 
for the temporal joint searcher and channel estimator the change in phase over time is 
the doppler frequency for the dimensional receptivity vector. The phase rotation speed, 

30 or frequency, of the dimensional receptivity vector, for the sampling window time 
instance can be interpreted as the doppler shift (DS). Thus, each dimensional 
receptivity vector corresponds to a separate doppler shift frequency. There are plural 
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possible frequencies for the dimensional receptivity vector, each of the plural possible 
frequencies corresponding to a different possible doppler shift for a wavefront. For the 
non-parametric techniques herein employed, the plural possible frequencies can be a 
continuous range of frequencies. For sake of differentiating the plural possible 
5 frequencies, the plural possible frequencies are each represented by a frequency index. 

[000138] For the temporal joint searcher and channel estimator, the channel 
estimate comprises, as mentioned before, a time of arrival (TOA) and doppler shift for 
each arriving wavefront in the sampling window (e.g., a channel coefficient mapped to 
a doppler shift). Therefore, the channel estimate may comprise a set (of one or more) 

10 pairs of data, each pair including a time of arrival (TOA) and a channel coefficient. 
The task for the temporal joint searcher and channel estimator is thus to locate a value 
or "tone" in antenna signal matrix 110 that best corresponds to an arriving wavefront, 
e.g., to locate a value or tone for each arriving wavefront in the sampling window. This 
task of locating a value or "tone" in antenna signal matrix 110 that best corresponds to 

15 an arriving wavefront can be accomplished by various techniques, including both 

parametric and non-parametric techniques. A Fast Fourier Transform (FFT) technique 
as discussed below is just one representative and illustrative example non-parametric 
type of correlator which can be utilized. 

[000139] Fig. 15 depicts example basic steps performed by an example correlator 
20 50-13A and correlator output analyzer 62-13A in conjunction with the Fast Fourier 
Transform (FFT) calculation. As step 15-1, the correlator 50-13A of Fig. 13A 
calculates Expression 5. 

Y(n,t) = FFT(n,X(n,t)) Expression 5 

wherein t is the sampling window time index; X(n,t) is the complex antenna matrix; 
25 and, n is the doppler frequency index. Each FFT calculation is thus a one dimensional 
FFT calculation on the baseband signal, and corresponds to a specific doppler shift 
frequency. 

[000140] The output of correlator 50-13A, i.e., the Y(n,t) values computed using 
Expression 1, are stored as correlator output values. The correlator output values can 
30 be stored, for example, in the correlator output value memory 52-13A of Fig. 13 A. 
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[000141] The correlator output analyzer 62-1 3 A of channel estimate (CE) generator 
60-13A searches the correlator output values and (as step 15-2) determines therefrom a 
maximum absolute value |Y(n,t)| max . This maximum absolute value |Y(n,t)| max is 
utilized by correlator output analyzer 62-13A to determine both the doppler shift (DS) 
5 and time of arrival (TO A) for an arriving wavefront. In particular, as step 15-3 the 
correlator output analyzer 62-13A chooses a sampling window time index t_max at 
which |Y(n,t)| max occurs to be the time of arrival of the arriving wavefront. In addition, 
as step 15-4 the correlator output analyzer 62-13A chooses the doppler index n_max at 
which |Y(n,t)| max occurs to determine the doppler shift (DS) of the arriving wavefront. 
10 An amplitude for the arriving wavefront is determined as correlator output analyzer 62- 
13A divides |Y(n,t)| max by the number of sets of pilot data comprising the antenna signal 
matrix (as step 15-5). 

[000142] Expression 5 and the steps of Fig. 15 represent a generic FFT calculation. 
In a CDMA-specific situation which utilizes a coding generator (such as coding 
15 generator 30 of Fig. 1), a comparable FFT calculation can be made using a refinement 
of Expression 5 such as that which appears as Expression 2, previously discussed, but 
applied for the temporal joint searcher and channel estimator rather than for the spatial 
joint searcher and channel estimator. 

[000143] As a result of operation of joint searcher and channel estimator 24- 13 A, 
20 an accurate channel estimate can be provided to the detector as a temporal signature. 
For each wavefront, the temporal signature includes the time of arrival (TOA) mapped 
to a doppler (frequency) shift. As explained below, the channel coefficient (CC) for 
each time of arrival and wavefront is derived from the doppler frequency shift. The 
time of arrival (TOA) and channel coefficient (CC) are applied to the detector as 
25 represented by lines 66-13A and 68- 13 A, respectively, in Fig. 13 A. 

[000144] As mentioned above, the channel coefficient (CC) for each wavefront is 
derived from the doppler frequency shift (DS). Recall that at step 15-4 the correlator 
output analyzer 62-2B chose the frequency index n_max at which |Y(n,t)| max occurs to 
represent the doppler shift frequency (DSF) of the arriving wavefront, with the chosen 
30 frequency index corresponding to a doppler shift (e.g., 0 ? , i.e., the derivative of 0). The 
channel impulse response vector (i.e., array propagation vector) x is therefore generated 
by detector interface 64-2B in accordance with Expression 6. 
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[000145] C[e> 2nn + H) , j 2ntT2 + H) , e^ 3 + H) , . . .J 2nm + H) ] Expression 6 

[000146] In Expression 6, C is the amplitude of the wavefront, f is the frequency of 
the signal (including doppler shift); T is the period time between two pilot 
symbols/sequences (which are assumed to be periodical, in analogy to the uniform 
5 array of the spatial embodiment), and H is a complex value of the signal at the first pilot 
symbol/sequence, H being arg(FFT max). For sake of simplicity, noise has been 
excluded from Expression 6, and C is assumed to be constant within the time TN. 

[000147] In the foregoing description it is the role of channel estimate (CE) 
generator 60-2A, and particularly detector interface 64-2A, to generate both a time of 

10 arrival (TOA) and a channel coefficient (CC), the channel coefficient being derived 

from the doppler shift, e.g., as above described in conjunction with Expression 6. In an 
alternate implementation of this and other embodiments described herein, the detector 
itself (such as detector 26 illustrated in Fig. 1), upon receiving the time of arrival 
(TOA) and doppler shift (DS) for each arriving wavefront, may have the intelligence to 

15 compute the channel coefficient for each wavefront from the corresponding direction of 
arrival (DO A) information. In such case, the time of arrival and direction of arrival are 
output by detector interface 64-13A to the detector. 

[000148] Thus, the joint searcher and channel estimator 24-13A looks at a discrete 
number of possible doppler frequency shifts, and picks the doppler frequency with the 
20 highest correlation (highest absolute value). 

[000149] Whereas the joint searcher and channel estimator of Fig. 13A includes a non- 
parametric correlator (e.g., a filter) which performs a Fast Fourier Transform (FFT) 
calculation, in other example embodiments the temporal joint searcher and channel 
estimator implements parametric techniques. As does the Fig. 13A embodiment, the 
25 spatial joint searcher and channel estimator 24-13B of Fig. 13B is shown along with its 
associated example antenna array 22-13B comprising an antenna element 22-13B-1 
which receives the successive sets of pilot data in the manner of Fig. 14. 

[000150] Similarly to the earlier described embodiment, joint searcher and channel 
estimator 24-13B can comprise an antenna signal matrix handling unit 40-13B, which 
30 in turn comprises antenna signal matrix generator 42-13B and antenna signal matrix 
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memory 44-13B, which function much in the manner previously described. For 
example, the complex baseband values stored in antenna signal matrix memory 44-13B 
can also be conceptualized as matrix 110, and as such has a sampling window time 
index. The antenna signal matrix 110 has been previously discussed in conjunction 
5 with Fig. 14, and is now also discussed with reference to Fig. 16A for sake of 
expounding the joint searcher and channel estimator 24-13B of Fig. 13B. 

[000151] The joint searcher and channel estimator 24-13B further comprises the 
parametric estimator 51-1 3B which outputs a parametric output estimation vector for 
storage in memory 52-13B. In addition, in similar manner as the preceding 
10 embodiment, joint searcher and channel estimator 24-13B comprises a channel estimate 
generator 60-13B which has parametric output estimation vector analyzer 62-13B and a 
demodulator interface 64-13B. Basic steps performed by parametric estimator 51-1 3B 
and parametric output estimation vector analyzer 62-13B of the joint searcher and 
channel estimator 24-1 3B of Fig. 13B are illustrated in Fig. 17. 

15 [000152] For each sampling window time index of the antenna signal matrix 110. As 
step 17-1, the parametric estimator 51-13B estimates, e.g., two parameters at each time 
instant: a temporal frequency parameter parameter and a temporal amplitude 
parameter. The temporal frequency parameter estimates the frequency the incident 
waves creates when arriving at the antenna for the consecutive pilot symbols. The 

20 temporal amplitude parameter estimates the amplitude of this frequency. The temporal 
frequency parameter and temporal amplitude parameter are considered to be a 
parameter pair and in Fig. 16B , they are illustrated as one parameter per sample along 
the sampling time index. 

[000153] As step 17-2 performed by joint searcher and channel estimator 24-1 3B, 
25 analyzer 62-13B finds certain "qualifying" values in parametric output estimation 

vector 120, i.e. maximum value of temporal amplitude vector. Each qualifying value of 
parametric output estimation vector 120 can correspond to an arriving wavefront for the 
sampling window. 

[000154] For each qualifying value, as step 17-3 the parametric output estimation 
30 vector analyzer 62-13B chooses a time of arrival (TO A) as corresponding to the 

sampling window time index t for the qualifying value, e.g., the time index at which the 
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maximum/qualifying absolute value of the parametric estimation output vector occurs 

[000155] Similarly, for each qualifying value, as step 17-4 the parametric output 
estimation vector analyzer 62-13B chooses a a doppler shift frequency (DS) as the 
temporal frequency parameter value at the time of arrival decided in 17-3. 

5 [000156] As step 17-5, the parametric estimation output vector analyzer 62-13B 

determines the amplitude as being the maximum/qualifying absolute value divided by 
the number of pilot data sets in the series. 

[000157] The joint searcher and channel estimator 24-13B thus looks for an 
optimum doppler (shift) frequency, and prepares a channel estimate which can be 

10 provided to the detector as a temporal signature. The temporal signature includes the 
time of arrival (TO A), as well as the doppler shift frequency (DSF) and amplitude. The 
channel coefficient (CC) for each time of arrival and wavefront is derived from the 
doppler shift (DS) in the manner described above with reference to Expression 6. The 
time of arrival (TO A) and channel coefficient (CC) are applied to the detector as 

15 represented by lines 66-13B and 68-13B, respectively, in Fig. 13B. 

[000158] It should be understood from the foregoing that information indicative of 
more than one incident wavefront may be seen in a sampling window. For example, 
with reference to the parametric output estimation vector 120 of Fig. 16B, the 
parametric output estimation vector analyzer 62-13B may see other high numbers and 

20 for each of those high numbers which qualify, an arriving wavefront may be 

ascertained. For example, if there were two high numbers, then the channel impulse 
response may reflect two arriving wavefronts. For each of the two arriving wavefronts 
the joint searcher and channel estimator would pick out both a time of arrival (TO A) 
and doppler shift frequency (DSF), as well as amplitude, which are mapped to two 

25 different channel coefficients, with these two different channel coefficients forming 
part of the channel estimate for the channel impulse response. 

[000159] The operation of the temporal searcher and channel estimator has been 
described above for one antenna element of an antenna array 22. It should be 
understood that the antenna array 22 may comprise plural antenna elements, and that 
30 the operations described above may be performed separately with respect to one or 



37 



more antenna elements of the array. Moreover, as described later, principles of the 
foregoing operation may be performed in a combined manner with respect to plural 
antennas of the antenna array. 

[000160] The temporal joint searcher and channel estimator and techniques of 
5 operation thereof as described above is particularly well suited for, but not limited to, a 
receiver unit which has only one antenna element, e.g., a mobile terminal with only one 
antenna. As indicated above, however, the temporal joint searcher and channel 
estimation techniques can be utilized separately but in parallel by plural antennas for a 
receiver. 

10 [000161] Consider, for example, the situation reflected in Fig. 1 1 in which the output of 
antenna element 22-13A-1 (or 22-13B-1) for a pilot data set Tl with the complex vector 
a M (and phase 0m); the output of the same antenna element for a pilot data set T2 with 
the complex vector a U2 (and phase 0i_ 2 ), and so forth. In this situation, the linear 
combinationof the complex antenna baseband signal and the doppler weight vectors Wi 

15 also has the effect of a summation, or coherent combination in the time domain, shown 
as summation function 100 in Fig. 12. By adding these complex vectors coherently, the 
temporal joint searcher and channel estimator increases the performance of the search 
and the channel estimate. 

[000162] In situations in which there is no doppler shift (e.g., the mobile terminal 
20 stands still or moves in a radial direction relative to the base station), the doppler shift 
frequency may be zero. In such cases the pilot data of the arriving wavefront(s) have 
essentially the same complex values. The situation of no doppler shift is just one 
special case of the generic operation of the temporal joint searcher and channel 
estimator described above. When the mobile starts to move, a doppler shift may occur, 
25 the temporal joint searcher and channel estimators obtains the doppler shift frequency, 
and thereby enhances the channel estimate. The channel estimate is enhanced by 
considering the doppler shift, regardless of the magnitude of the doppler shift. 

[000163] The non-parametric FFT-type correlator and the parametric estimator 
techniques illustrated above, e.g., by Fig. 13A and Fig. 13B, respectively, are only two 
30 example techniques for finding the values or 11 tones' 1 in antenna signal matrix 110. 
Other parametric approaches are described by or understood from Stocia, Petre and 



38 



Moses, Randolph, Introduction To Spectral Analysis, ISBN-01 3-258419-0, Prentice 
Hall, which is incorporated by reference in its entirety, particularly Chapter 4 thereof. 

[000164] SPATIAL-TEMPORAL JOINT SEARCHERS/ESTIMATORS 

[000165] In some further embodiments, which combine features from both the 
5 spatial and temporal embodiments discussed above, plural antenna elements of an 

antenna array provide respective plural series of signals for successive sets of pilot data. 
The joint searcher and channel estimators of these further embodiments essentially 
concurrently consider the plural series of signals provided by the plural antennas for 
determining both a time of arrival and channel coefficient. 

10 [000166] By concurrently considering the signals provided by plural antennas, the 
channel estimate takes into consideration direction of arrival in determining the time of 
arrival and channel coefficient. By concurrently considering the series of signals 
provided by each antenna, in which each series comprises successive sets of pilot data, 
the channel estimate further takes into consideration a frequency shift which may be a 

15 doppler shift (occasioned by relative movement of a transmitter and the receiver or of 
an object in a field between the transmitter and receiver). The channel estimate is 
performed by considering spatial and temporal domain jointly and concurrently. 

[000167] Since it processes the series of signals from plural antennas, with each 
series comprising successive sets of pilot data, the joint searcher and channel estimator 

20 is considered a three dimensional unit. A first dimension is with reference to a time 
index of a sampling window, i.e., a sampling window time index. A second dimension 
is a spatial dimension imparted by the spacing of the plural antennas of the array. This 
spatial dimension, which involves essentially simultaneous and concurrent processing 
together of signals from the plural antennas for the array in order to determine the time 

25 of arrival and channel coefficient, bestows on the joint searcher and channel estimator 
the distinction of being a "spatial" joint searcher and channel estimator. A third 
dimension is a temporal dimension imparted by the time interval reflected by the 
successive sets of pilot data. This temporal dimension, which involves essentially 
simultaneous and concurrent processing together of signals for each of the successive 

30 sets of pilot data in order to determine the time of arrival and channel coefficient, 
bestows on the joint searcher and channel estimator the distinction of being a 
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"temporal" joint searcher and channel estimator. In view of being both a spatial and 
temporal joint searcher and channel estimator, the joint searcher and channel estimator 
is also referred to as a "combined" spatial/temporal joint searcher and channel 
estimator, or spatio/temporal joint searcher and channel estimator. 

5 [000168] Concurrent consideration of the plural series of signals may be either in a 
three dimensional essentially concurrent mode or a sequenced mode. The three 
dimensional essentially concurrent mode involves a single step determination of the 
time of arrival and channel coefficient by simultaneously considering signals from all 
antennas of the array for all of the plural series. The sequenced mode involves a two 

10 step determination of the time of arrival and channel coefficient. In the sequenced 
mode, a first step comprises determining a time of arrival and direction of arrival by 
concurrently considering the plural signals provided by the plural antennas for a first of 
the plural series. A second step of the sequenced mode comprises refining the 
estimation of the channel coefficient based on doppler shift by concurrently considering 

15 the elements of the plural series having the direction of arrival determined in the first 
step. This procedure could also be performed the other way around: first determining 
the time of arrival and Doppler shift and then refining the channel estimate by 
concurrently considering the elements of the plural series having the Doppler shift 
determined in the first step. 

20 [000169] Fig. 18A illustrates one example embodiment of a spatio-temporal joint 
searcher and channel estimator 24- 13 A, as well as an associated example antenna array 
22-1 8A. The antenna array 22-1 8A includes, by way of non-limiting example, four 
antenna elements 22-18A-1 through 22-18A-4. While the antenna elements 22-18A-1 
through 22-1 8 A-4 are shown as forming a uniform linear array (ULA), it should be 

25 understood that antenna configurations other than a uniform linear are possible, and that 
the number of antenna elements in the antenna array may vary (e.g., the number of 
antenna elements is not limited to four). After suitable radio frequency processing, the 
signals obtained from the antenna elements are each applied as baseband signals to joint 
searcher and channel estimator 24- 18 A, as well as to a detector (not illustrated in Fig. 

30 18 A). 

[000170] The joint searcher and channel estimator 24-18A comprises an antenna signal 
matrix handling unit 40-18A. In one particular example manifestation, antenna signal 
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matrix handling unit 40-18A includes antenna signal matrix generator 42-18A and 
antenna signal matrix memory 44-18A. A matrix analyzer, which for the non- 
parametric technique of Fig. 18A can be a correlator 50- 18 A, operates on complex 
values stored in antenna signal matrix memory 44- 18 A. The correlator 50-1 8 A 
5 preferably comprises a filter. The correlator 50-1 8 A generates certain output values, 
which may be stored, e.g., in correlator output value memory 52- 18 A. The joint 
searcher and channel estimator 24-18A further comprises a channel estimate (CE) 
generator 60-18A. In the illustrated example embodiment, the channel estimate (CE) 
generator 60-18A comprises a correlator output analyzer 62-18A and a detector 
10 interface 64- 18 A. The detector interface 64-18A generates, for each wavefront, a 

channel estimate which includes both a time of arrival (TO A) and a channel coefficient 
(CC). In Fig. 18 A, the time of arrival and channel coefficient output by detector 
interface 64 are applied to the detector on lines 66-1 8 A and 68- 18 A, respectively. 

[000171] In the joint searcher and channel estimator 24-18A of Fig. 18A, for each series 
15 of sets of pilot data (represented by pilot data sets Tl - T4), the antenna matrix handling 
unit 40-18A samples the signals from each antenna element. Using the sampled 
signals, antenna signal matrix generator 42-18A generates an antenna signal matrix 
such as antenna signal matrix 130 illustrated in Fig. 19. The antenna signal matrix 130 
may be stored in any convenient fashion, such as antenna matrix memory 44- 18 A. 

20 [000172] The antenna signal matrix 130 is a three dimensional functionally dependent 
matrix. In other words, complex samples are stored in antenna signal matrix 130 as a 
function of three different indexes. For the antenna signal matrix 130 shown in Fig. 19, 
a first index is a sampling window time index, illustrated along the X axis of Fig. 19. 
For embodiments which utilize spreading codes or similar codes, the first index may be, 

25 for example, a chip index. Thus, the sampling window time index points to a time in 
the sampling window relative to a start of the sampling window. 

[000173] In the antenna signal matrix 130 of Fig. 19, a second index, shown along the 
Y axis, is an antenna index. The antenna index points to a different row of the antenna 
signal matrix 130, each row being associated with a different antenna element in 
30 antenna array 22. Fig. 19 shows four rows in antenna signal matrix 130 for consistency 
with the previous examples of an antenna array comprising four antenna elements. It is 
reiterated, however, that the number of antennas in an antenna array, and thus the 
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number of rows in antenna signal matrix 130 and the maximum value of the antenna 
index, can vary from receiver to receiver, and that the choice of four antenna is only 
illustrative for sake of example. 

[000174] In the antenna signal matrix 130 of Fig. 19, a third index, shown along the Z 
5 axis, is a pilot set index. The pilot set index indicates which one of the sets of pilot data 
the sample was obtained. In other words, a pilot set index = Tl indicates that the 
sample was obtained from pilot set Tl; a pilot set index = T2 indicates that the sample 
was obtained from pilot set T2; and so on as depicted by the arrows which connect the 
matrix 110 with the received signal depiction with its illustrative successive sets of pilot 
10 data. As can be seen, the pilot set index points to a different plane of the antenna signal 
matrix 110, each plane being associated with a different set of pilot data. 

[000175] Fig. 19 shows four planes in antenna signal matrix 130 for consistency with 
the illustrated example wherein the antenna signal matrix encompasses four successive 
sets of pilot data. The number of sets of pilot data subsumed in a given antenna signal 

15 matrix, and thus the maximum value of the pilot set index, can vary from receiver to 
receiver, so that the present example's choice of four sets of pilot data is only 
illustrative for sake of example. In general, the choice of the number of sets of pilot 
data to be apprehended simultaneously by a spatio-temporal spatial/temporal joint 
searcher and channel estimator depends on how quickly the doppler is expected to 

20 change. The number of taps/incident waves depends on the multipath. In other words 
in an open space we have one direct path and thus only one channel/tap coefficient in 
the channel impulse response. 

[000176] For sake of simplicity, the complex values stored in antenna signal matrix 
130, including the complex values obtained from the antennas, are not illustrated in Fig. 
25 19. Such complex values would be illustrated in a fourth dimension. 

[000177] In conjunction with the antenna signal matrix 130 of Fig. 19, and particularly 
a WCDMA case in which spacing of the antenna elements in the antenna array is not 
too far apart, the plane wavefront arriving at the antenna array can be considered to 
arrive in the same sampling window time index (or chip index). 
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[000178] Assuming that the wavefront arrives at the antenna elements at different 
times (the time differences are small in comparison to the sampling time interval), the 
complex values stored for each column of the antenna signal matrix 130 of Fig. 19 have 
differing phase (e.g., 9) values in each row of the column. For a uniformly spaced array 
5 of antenna elements, the phase difference is essentially the same between adjacent rows 
of the same column (although noise may be a factor). But whatever the spacing, the 
rate of change of the phase with respect to time (time of travel of the approaching 
wavefront) is the phase rotation speed, or frequency, for the vector formed by the 
column, as previously explained. This per column frequency can be interpreted as a 
10 direction of arrival (DO A). There are plural possible frequencies for the columns of 

antenna signal matrix 130, with each of the plural possible frequencies corresponding to 
a possible direction of arrival (DOA) of a wavefront. The plural possible direction of 
arrival frequencies are represented by a frequency index "ni". 

[000179] In a similar manner, for each slice of antenna signal matrix 130 along the 
15 Z" direction the complex values have differing phase (e.g., 0) values. The Z-aligned 
elements of differing "Z" planes of antenna signal matrix 130 have differing phase 
values in view of a possible doppler shift as detected by the different sets of pilot data 
as gathered over plural sets of pilot data in a series. The rate of change over time of the 
phase along the Z direction between successive sets of pilot data is a frequency which is 
20 associated with the doppler shift. There are plural possible frequencies for the Z slices 
of antenna signal matrix 130, with each of the plural possible frequencies 
corresponding to a possible doppler shift (DS) for wavefront. The plural possible 
doppler shift frequencies are represented by a frequency index "n 2 ". 

[000180] The channel estimate generator 60-18A (see Fig. 18 A) seeks to develop a 
25 "composite" channel estimate based on the complex values stored in antenna signal 

matrix 130. As mentioned before, since an antenna array such as antenna array 22-18A 
has plural antenna elements, there are a corresponding plurality of channels through 
which wavefronts are received, and accordingly there could also be a separate channel 
impulse response or separate channel estimate for each of the plural channels. But by 
30 storing the complex samples in antenna signal matrix 130 in the manner aforedescribed, 
and by concurrently finding the time of arrival (TOA) and channel coefficients over the 
entire antenna signal matrix 130, the channel estimate generator 60-1 8 A provides a 
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channel estimate which encompasses all channels for all antenna elements and for this 
reason is known as a "composite" channel estimate. 



[000181] The composite channel estimate comprises, as mentioned before, a time 
of arrival (TO A) and channel coefficient for each arriving wavefront in the sampling 

5 window (e.g., a channel coefficient mapped to a time of arrival (TOA)). Therefore, the 
channel estimate may comprise a set (of one or more) pairs of data, each pair including 
a time of arrival (TOA) and channel coefficient. The task for correlator 50-1 8 A is thus 
to locate a value or "tone" in antenna signal matrix 130 that best corresponds to an 
arriving wavefront, e.g., to locate a value or tone for each arriving wavefront in the 

10 sampling window. 

[000182] The task of locating, in an antenna signal matrix such as antenna signal 
matrix 130, a value or "tone" that best corresponds to an arriving wavefront can be 
accomplished by various techniques, including both parametric and non-parametric 
techniques. A Fast Fourier Transform (FFT)) performed in a three dimensional 
15 essentially concurrent mode is discussed below in conjunction with as just one 
representative and illustrative example of a non-parametric technique wherein 
correlator 50-18A is utilized. 

[000183] Fig. 20 depicts example basic steps performed by an example correlator 
50-1 8 A and analyzer 62-18A in conjunction with the Fast Fourier Transform (FFT) 
20 calculation. In conjunction with Fig. 20, Fig. 21 shows an antenna signal matrix; a 
doppler weight and antenna weight vector; and a non-parametric estimation output 
vector for an example embodiment of a spatio-temporal joint searcher and channel 
estimator which operates in a three dimensional essentially concurrent mode. As step 
20-1, the correlator 50-18A of Fig. 18A calculates Expression 8. 

25 Y(n l9 n 2 ,t) = FFT(n 1 ,n 2 ,X(:,:t)) Expression 8 

In Expression 8, t is the sampling window time index; X(:,:,t) is the complex antenna 
matrix (with the colon representing all antenna indexes for one sampling window 
time index); n x is the direction of arrival frequency index; and n 2 is the doppler shift 
index. Each FFT calculation is thus a two dimensional FFT calculation on the 
30 baseband signal, corresponding both to a specific direction of arrival (as depicted by the 
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frequency index n^ and to a specific doppler shift (as depicted by the frequency index 
n 2 ). 

[000184] The output of correlator 50-18A, i.e., the Y(nj,n 2 ,t) values computed 
using Expression 8, are stored as correlator output values. The correlator output values 
5 can be stored, for example, in the correlator output value memory 52-18A of Fig. 18A. 

[000185] The correlator output analyzer 62-18A of channel estimate (CE) generator 
60-18A searches the correlator output values Y(ni,n 2 ,t) and (as step 20-2) determines 
therefrom a maximum absolute value |Y(n 1 ,n 2 ,t)| max . This maximum absolute value 
|Y(ni,n 2 ,t)| max is utilized by correlator output analyzer 62-18A to determine both the 

10 direction of arrival (DO A) and time of arrival (TO A) for an arriving wavefront seen in 
the sampling window. In particular, as step 20-3 the correlator output analyzer 62-18A 
chooses a sampling window time index t_max at which |Y(n l5 n 2 ,t)| max occurs to be the 
time of arrival of the arriving wavefront. In addition, as step 20-4 the correlator output 
analyzer 62-1 8 A chooses the frequency index n t _max at which |Y(ni,n 2 ,t)| max occurs to 

15 determine the direction of arrival (DOA) of the arriving wavefront. Further, as step 20- 
5 the correlator output analyzer 62-18A chooses the index n^max at which 
|Y(ni,n 2 ,t)| max occurs to determine the doppler shift of the arriving wavefront. An 
amplitude for the arriving wavefront is determined as the correlator output analyzer 62- 
18A divides |Y(ni,n 2 ,t)| max by the product of the number of antennas comprising the 

20 antenna array and the number of sets of pilot data included in the matrix 130 (as step 
20-6). 

[000186] Expression 8 and the steps of Fig. 20 represent a generic FFT calculation. 
In a CDMA-specific situation which utilizes a coding generator (such as coding 
generator 30 of Fig. 1), a comparable FFT calculation can be made using a refinement 
25 of Expression 8 which appears as Expression 9. 

Y(nj,n 2 ,t) = 2 q *FFT(n 1 ,n 2 ,X(:,:,t)), j = 1,K Expression 9 

Expression 9 is understood from Expression 1, it being further mentioned that Cj is a 
coding sequence symbol value j; and K is a length of the coding sequence. 



[000187] As a result of operation of spatio-temporal joint searcher and channel 
estimator 24- 18 A, an accurate channel estimate can be provided to the detector as a 
spatio-temporal spatial and temporal signature. The spatial signature includes the 
direction of arrival; the temporal signature includes the doppler shift. The channel 
5 coefficient (CC) for each time of arrival and antenna element is derived from the 

direction of arrival (DOA)and the doppler shift. The time of arrival (TOA) and channel 
coefficient (CC) are applied to the detector as represented by lines 66-18A and 68- 18 A, 
respectively, in Fig. 18 A. 

[000188] As mentioned above, the channel coefficient (CC) for each wavefront is 
10 derived from the direction of arrival (DOA) and doppler shift (DS). Recall that at step 
18-4 analyzer 62- 18 A chose the frequency index n u max at which |Y(ni,n 2 ,t)| max occurs 
to represent the direction of arrival (DOA) of the arriving wavefront, with the chosen 
frequency index corresponding to a direction of arrival (e.g., 0). Further, analyzer 62- 
18A chose the frequency index n 2 _max at which |Y(ni,n 2 ,t)| max occurs to represent the 
15 doppler shift of the arriving wavefront, with the chosen frequency index corresponding 
to a doppler shift. The channel impulse response vector (i.e., array propagation vector) 
x is therefore generated by detector interface 64-18A in accordance with Expression 10 
(for identical isotropic antenna elements). 



X =[(l,e (jkd * Sin0) , e 01 " 1 * 2 ^, . e 0kd*(K-l)sin9)j* CO; 
20 (l 5e (ikd*sin8) e Gkd*2sin8) e Gkd*(K-l)sin8)] +cl . 



(l,e« kd * sin9 >, e akd * 2sin0) , . . . e » 1)sin0) ]*CN Expression 10 

In Expression 10, CN = e i 2nnN = H ^ with H and other parameters being as previously 
defined. 

[000189] In the foregoing description it is the role of channel estimate (CE) 
25 generator 60- 18 A, and particularly detector interface 64- 18 A, to generate both a time of 
arrival (TOA) and a channel coefficient (CC), the channel coefficient being derived 
from the direction of arrival and doppler shift, e.g., as above described in conjunction 
with Expression 11. In an alternate implementation of this and other embodiments 
described herein, the detector itself (such as detector 26 illustrated in Fig. 1), upon 
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receiving the time of arrival (TO A), the direction of arrival (DO A), and the doppler 
shift for each arriving wavefront, may have the intelligence to compute the channel 
coefficient for each wavefront from the corresponding direction of arrival (DOA) and 
doppler shift information. In such case, the time of arrival, direction of arrival, and 
5 doppler shift are output by detector interface 64 to the detector. 

[000190] The operation of the correlator 50-18A in calculating Expression 8 or 
Expression 9 is an example of a three dimensional essentially concurrent mode, since 
evaluation of Expression 8 (or Expression 9 for a WCDMA implementation) involves a 
single step determination of the time of arrival and channel coefficient by simultaneous 
10 considering signals from all antennas of the array for all of the plural series. In other 
words, in the illustrated example of the three dimensionally essentially concurrent 
mode, the Fast Fourier Transform (FFT) of Expression 8 or Expression 9 had three 
arguments: n l5 n 2 , and X(:,:t), so that the FFT operated on all arguments essentially 
simultaneously. 

15 [000191] In contrast to the three dimensional essentially concurrent mode, the 
sequenced mode involves a two step determination of the time of arrival and channel 
coefficient. In a first alternative way for implementing the sequenced mode, a first step 
comprises determining a time of arrival and direction of arrival by concurrently 
considering the plural signals provided by the plural antennas for a first of the plural 

20 series. For example, the first step of the first alternative of the sequenced mode can 
involve calculating a FFT such as that of Expression 1 (or, for WCDMA, Expression 
2). From the results of the first step or first FFT calculation, a time of arrival (TO A) 
and tentative channel coefficient are determined. Then, as a second step of the first 
alternative of the sequenced mode, the tentative channel coefficient is refined by taking 

25 into consideration a possible frequency shift (e.g., doppler shift) by further considering 
the elements of the plural series having the direction of arrival determined in the first 
step. In a second alternative way of implementing the sequenced mode, the order of the 
steps is essentially reversed: first the FFT is performed in the temporal domain to 
decide time of arrival and tentative channel coefficient; and secondly the tentative 

30 channel coefficient is refined by FFT in the spatial domain. 



[000192] The procedures of the first alternative implementation of the sequenced 
mode for the non-parametric technique are illustrated in Fig. 22 A and Fig. 22B in 
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conjunction with Fig. 23. Fig. 22A and Fig. 22B as a diagrammatic view of an antenna 
signal matrix; an antenna weight vector; and a non-parametric estimation output vector 
for an example embodiment of a sequential spatio-temporal joint searcher and channel 
estimator. In Fig. 22A, the FFT operates on the spatial domain and calculates the FFT 
5 (illustrated by the FFT vector Wi) for the antenna matrix for each time interval. The 
time of arrival is choosen by picking the direction of arrival index and time index with 
the highest absolute value. If this index does not coincide for all time intervals, the 
index can be chosen with some method, e.g. majority decision. 

[000193] After having chosen the time of arrival index and direction of arrival 
10 index, these FFT -processed samples are further FFT-processed by a FFT calculation in 
the temporal domain (illustrated by FFT frequency vector Wj). Fig. 22B shows the 
spatially filtered samples for the identified time of arrival and direction (marked as grey 
in the figure) are filtered with the temporal vectors. After the second FFT processing, 
the channel estimate is created from the sample with the highest magnitude. Step 23-1 
15 through step 23-7 of Fig. 23 also describe the procedure of the first alternative 
implementation of the sequenced mode. 

[000194] The procedures of the second alternative implementation of the sequenced 
mode for the non-parametric technique are illustrated in Fig. 24A and Fig. 24B in 
conjunction with Fig. 25. Fig. 24A and Fig. 24 B show an antenna signal matrix; a 

20 Doppler weight vector; and a non-parametric estimation output vector. In Fig. 24A the 
FFT operates on the temporal domain and calculates the FFT (illustrated by the FFT 
vector Wj) for the antenna matrix for each time interval. The time of arrival is chosen 
by picking the doppler index of and time index with the highest absolute value. If this 
index does not coincide for all time intervals, the index can be chosen with some 

25 method, e.g. majority decision. After having chosen the time of arrival index and 
doppler index, these FFT-processed samples are further FFT-processed by a FFT 
calculation in the spatial domain (illustrated by FFT frequency vector Wi). Fig. 24B 
shows the spatially filtered samples for the identified time of arrival and Doppler shift 
(marked as grey in the figure) are filtered with the spatial vectors. After the second 

30 FFT processing, the channel estimate is created from the sample with the highest 

magnitude. Step 25-1 through step 25-7 of Fig. 25 also describe the procedure of the 
second alternative implementation of the sequenced mode. 
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[000195] Whereas the joint searcher and channel estimator of Fig. 18A includes a non- 
parametric type correlator (e.g., a filter which performs a Fast Fourier Transform (FFT) 
calculation), in other example embodiments the joint searcher and channel estimator 
implements parametric techniques. As does the Fig. 18A embodiment, the parametric 
5 temporal joint searcher and channel estimator 24-18B of Fig. 18B is shown along with 
its associated example antenna array 22-18B. Again by way of example, antenna array 
22-18B includes four antenna elements 22-18B-1 through 22-18B-4. The signals 
obtained from the antenna elements are each applied to joint searcher and channel 
estimator 24-18B, as well as to a detector (not illustrated in Fig. 18B). 

10 [000196] Similarly to the earlier described embodiment, joint searcher and channel 
estimator 24-18B can comprise an antenna signal matrix handling unit 40-18B, which 
in turn comprises antenna signal matrix generator 42-18B and antenna signal matrix 
memory 44-18B, which function much in the manner previously described. For 
example, the complex baseband values stored in antenna signal matrix memory 44-18B 

15 can also be conceptualized as matrix 130, and as such has a sampling window time 

index. The antenna signal matrix 80 has been previously discussed in conjunction with 
Fig. 19. 

[000197] The joint searcher and channel estimator 24-1 8B further comprises a 
parametric estimator 51-1 8B which produces a parametric estimation output vector. In 
20 addition, in similar manner as the preceding embodiment, joint searcher and channel 
estimator 24-18B comprises a channel estimate generator 60-18B which has parametric 
output estimation vector analyzer 62-18B and a demodulator interface 64-18B 

[000198] Fig. 26 shows an antenna signal matrix and a parametric estimation output 
vector for an example embodiment of a spatio-temporal joint searcher and channel 
25 estimator. As with the non-parametric techniques, the parametric techniques can be 
implemented either in a three dimensional essentially concurrent mode or in a 
sequenced mode, with the sequence mode having two alternative implementations. 

[000199] Fig. 27 shows basic, representative steps involved in a parametric three 
dimensional essentially concurrent mode. Step 27-1 shows the joint searcher and 
30 channel estimator 24-18B producing a parametric estimation output vector. Then, as 
step 27-2, analyzer 62-18B finds the "qualifying" values in the parametric estimation 



output vector. 



[000200] For each qualifying value, as step 27-3 the parametric output estimation 
vector analyzer 62-18B chooses a time of arrival (TOA) as corresponding to the 
sampling window time index t for the qualifying value, e.g., the time index at which the 
5 maximum/qualifying absolute value of the parametric estimation output vector occurs. 

[000201] For each qualifying value, as step 27-4 the parametric output estimation 
vector analyzer 62-18B chooses a spatio-temporal frequency parameter corresponding 
to the spatio-temporal frequency for the maximum/qualifying absolute value of the 
parametric estimation output vector. 

10 [000202] As step 27-5, the parametric estimation output vector analyzer 62-13B 

determines the amplitude as the value for the spatio-temporal amplitude value for the 
time of arrival decided in step 27-2. 

[000203] It should be understood from the foregoing that information indicative of 
more than one incident wavefront may be seen in a sampling window. For example, 
15 with reference to the parametric estimate output vector 140 of Fig. 26, the parametric 
output estimation vector analyzer 62-18B may see other (e.g., plural) high numbers, and 
for each of those high numbers which qualify, an arriving wavefront may be 
ascertained. 

[000204] The procedures of the first alternative implementation of the sequenced 
20 mode for the parametric technique are illustrated in Fig. 28A and Fig. 28B in 

conjunction with Fig. 29. Fig. 28A and Fig. 28B depict a parametric, sequential spatio- 
temporal joint searcher and channel estimator for this first alternative implementation. 
In Fig. 28A and Fig. 28B the parametric approach first operates on the spatial domain 
and calculates the spatial frequency parameters for each time instant over the time 
25 transmission intervals. The time of arrival is chosen by picking the spatial frequency 
amplitude value with the highest absolute value. The direction of arrival, DOA, is the 
value of the spatial frequency parameter. If this time of arrival does not coincide for all 
the time intervals, the time of arrival can be chosen with some method, e.g. majority 
decision. As shown in Fig. 28B, after having chosen the time of arrival index and 
30 direction of arrival, these samples are processed by the parametric approach applied in 
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the temporal domain. After the second processing, the channel estimate is created from 
the temporal parameters. Step 29-1 through step 29-5 of Fig. 29 also describe the 
procedure of the first alternative implementation of the parametric sequenced mode. 

[000205] Fig. 30A and Fig. 30B show a parametric, sequential spatio-temporal 
joint searcher and channel estimator for a second alternative implementation of a 
parametric, sequential spatio-temporal joint searcher and channel estimator. In Fig. 
30A and Fig. 30B the paremetric approach first operates on the temporal domain and 
calculates the temporal frequency parameters for each time instant over the time 
transmission intervals. The time of arrival is chosen by picking the temporal frequency 
amplitude value with the highest absolute value. The Doppler shift frequency, DSF is 
the value of the temporal frequency parameter. If this time of arrival does not coincide 
for all the time intervals, the time of arrival can be chosen with some method, e.g. 
majority decision. As shown in Fig. 30B, after have chosen the time of arrival index 
and DSF, these samples are processed by parametric approach applied in the spatial 
domain. After the second processing, the channel estimate is created from the spatial 
parameters. Step 31-1 through step 31-7 of Fig. 31 also describe the procedure of the 
second alternative implementation of the parametric sequenced mode. 

[000206] The non-parametric FFT-type correlator and the parametric linear 
combination logic techniques illustrated above are only two example techniques for 
finding the values or "tones" in antenna signal matrix 130 which are associated with 
arriving wavefronts. Other parametric approaches are described by or understood from 
Stocia, Petre and Moses, Randolph, Introduction To Spectral Analysis, ISBN-013- 
258419-0, Prentice Hall, which is incorporated by reference in its entirety, particularly 
Chapter 4 thereof. 

[000207] The spatio-temporal joint searcher and channel estimator and techniques 
of operation thereof as described above are suitable for any receiver unit which has 
plural receiving antennas. Thus, the spatial joint searcher and channel estimator is 
particularly well suited for, but not limited to, a base station which has plural antennas. 
The spatio-temporal joint searcher and channel estimator and techniques of operation 
thereof also encompasses mobile terminals which have plural antennas. 
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[000208] The joint searcher and channel estimators thus employ a multi- 
dimensional and optimum detection and estimation approach. The multi-dimensional 
joint searcher and channel estimators typified by those described herein have better 
performance than traditional one dimensional searchers. The multi-dimensional joint 
searcher and channel estimators have a greater SNIR for detecting time of arrival, 
which increases the probability of that the correct time of arrival will be ascertained. 
This, in turn, leads to a better channel estimate. 

[000209] In terms of implementation, the blocks, units, and functionalities of the 
differing embodiments of the joint searcher and channel estimator herein described can 
various forms. For example, those skilled in the art will appreciate that one or more of 
the functionalities of the joint searcher and channel estimator can be implemented using 
individual hardware circuits, using software functioning in conjunction with a suitably 
programmed digital microprocessor or general purpose computer, using an application 
specific integrated circuit (ASIC), and/or using one or more digital signal processors 
(DSPs). Furthermore, the functionalities of the joint searcher and channel estimator 
need not be delineated specifically in the manners illustrated, it being understood (for 
example) that the functionalities can be distributed, combined, subdivided, or otherwise 
rearranged for accomplishing essentially the same results. 

[000210] Use and operation of the joint searcher and channel estimators is not 
confined to WCDMA transmission, although in some instances WCDMA has been 
described above as an example environment of implementation. The principles, 
techniques, methods, and apparatus described herein can be adapted or augmented for 
compatibility with various types of networks, not only WCDMA, but other networks as 
well (such as GSM, for example). 

[00021 1] In the foregoing, it will be appreciated that other aspects of wireless 
receiver structure and operation which are tangential to matters described above have 
been omitted for clarity. Such aspects, well understood by persons skilled in the art, 
include without limitation pulse shaping, sampling frequency, time jitter, time 
alignment, demodulation, inter symbol interference (ISI), and co channel interference 
(CCI). 
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[000212] While the invention has been described in connection with what is 
presently considered to be the most practical and preferred embodiment, it is to be 
understood that the invention is not to be limited to the disclosed embodiment, but on 
the contrary, is intended to cover various modifications and equivalent arrangements 
5 included within the spirit and scope of the appended claims. 



